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1.    Reasons  for  a  Policy  Model 

Much  has  been  made  of  an  energy  "crisis"  in  the  Urited  States,  and  the 
need  for  planning  to  meet  "needs"  of  consumers  in  the  coming  decades. 
Economists  have  not  joined  the  legions  of  those  proclaiming  "emergency," 
except  perhaps  in  the  case  of  natural  gas.   In  this  industry,  there  is  a 
"crisis"  in  the  economist's  sense  of  non-clearing  of  demands.   There  exists 
a  shortage  of  natural  gas  in  this  country,  to  the  extent  that  non-delivery 
on  contract  sales  was  approximately  3  percent  of  total  sales  of  the  inter- 
state pipelines  in  1971-1972,  and  2  percent  of  national  sales  by  all  sources 
[4,  Chapter  4].   This  crisis  has  been  precipitated  by  direct  control  of  gas 
field  contract  prices  by  the  Federal  Power  Commission  in  the  1960's  [18]. 
Policy  changes  are  needed  to  solve  this  "crisis." 

One  of  the  purposes  of  a  policy  model  of  natural  gas  will  be  to  test 
that  assertion.   By  building  an  econometric  model  of  gas  markets  with  expli- 
cit policy  controls,  and  simulating  over  time,  we  expect  to  be  able  to 
examine  the  effects  of  past  regulatory  policy,  as  well  as  to  perform  ex- 
periments to  find  the  effects  of  alternative  regulatory  policies.   Finally, 
by  performing  simulations  of  the  model  into  the  future  using  different 
assumptions  about  what  regulatory  practices  of  the  Federal  Power  Commission 
might  be,  we  plan  to  derive  a  set  of  alternative  forecasts  for  excess  demand 
at  both  the  producer  and  the  wholesale  levels. 

Most  previous  econometric  studies  of  natural  gas  have  investigated 
either  demand  or  supply  of  gas,  but  have  neglected  the  simultaneous  inter- 
action of  these  two  sides  of  markets.   Balestra  [2],  for  example,  in  his 


classic  study  of  the  demand  for  natural  gas  by  residential  and  commercial 
consumers,  assumed  a  perfectly  elastic  supply  curve  for  production.   This 
assumption  was  probably  justified  during  the  1950's  and  1960's  since  pro- 
duction of  gas  for  final  consumers  took  place  on  an  "as  needed"  basis  from 
large  stocks  of  reserves.   This  modeling  approach  would  not  continue  to  be 
valid  during  the  1970' s,  however,  as  total  gas  demand  exceeds  the  constraints 
on  production  imposed  by  smaller  reserve  levels.   The  supply  studies  of 
Erickson  and  Spann  [6]  and  Khazzoom  [14],  similarly,  are  certainly  admirable 
attempts  at  defining  and  testing  some  of  the  relationships  that  exist  in 
the  gas  industry,  particularly  those  accounting  for  reduced  reserve  levels 
under  price  controls.   But,  to  the  extent  that  policies  are  changed  in  the 
future  so  that  markets  clear,  and  demand  is  once  again  observed,  models  of 
only  the  supply  side  of  markets  will  be  inadequate  to  represent  the  effects 
of  policy.   If  the  "crisis"  is  to  be  ameliorated,  in  other  words,  the  pro- 
duction and  reserve  supply  levels  of  the  industry  have  to  be  analyzed  as  a 
simultaneous  system. 

The  starting  point  for  this  model  is  a  recent  simultaneous  equation 
study  by  MacAvoy  [18]  which  focuses  primarily  on  demonstrating  how  price 
regulation  at  the  wellhead  has  led  to  reserve  shortages  and  a  consequent 
inability  of  field  reserve  markets  to  be  cleared  of  excess  demand.   His 
model  treats  only  the  market  for  new  reserves  (i.e.,  at  the  wellhead)  with 
the  pre-regulation  wellhead  price  acting  as  a  demand  variable  that  is  deter- 
mined in  part  by  exogenous  variables  (population,  income,  all-fuels  price 
index)  that  in  fact  determine  the  demand  for  ultimate  gas  consumption.   It 


serves  to  demonstrate  that  the  production  market,  and  excess  demand  in 
that  market,  can  be  at  least  roughly  described  by  a  simple  simultaneous 
equation  econometric  model. 

The  policy  model  developed  here  treats  simultaneously  the  regulation- 
induced  disequilibrium  in  the  field  market  for  reserves  (gas  producers 
selling  to  pipeline  companies  at  the  wellhead  price)  and  equilibrium  in 
the  wholesale  market  for  production  (pipeline  companies  selling  to  public 
utilities  and  to  industrial  consumers  at  different  wholesale  prices) .   The 
linking  of  these  two  markets  is  an  important  characteristic  of  the  natural 
gas  industry:   production  in  the  wholesale  market  is  a  determinant  of  pipe- 
lines' demand  for  gas  reserves  in  the  field  market,  and  price  on  new  reserve 
contracts  in  the  field  market  is  a  determinant  of  pipeline  delivery  costs 
and  thus  price  in  the  wholesale  market.   Disequilibrium  in  the  wholesale 
market  occurs  when  meeting  total  consumers  demands  would  result  in  total 
production  that  would  exceed  the  constraints  imposed  by  reserve  levels  in 
the  field  market. 

The  structure  of  the  model  is  described  in  some  detail  in  the  next 
section,  as  a  set  of  econometric  relationships  between  several  aggregated 
market  and  policy  variables.   Variables  which  are  endogenous  to  the  field 
market  include,  on  the  supply  side,  non-associated  and  (oil)  associated 
discoveries  of  reserves,  extensions  and  revisions  of  reserves,  and  wells 
drilled.   On  the  demand  side,  wellhead  prices  would  be  endogenous  if  demands 
were  cleared,  but  after  ceiling  prices  were  set  by  the  F.P.C.  in  the  1960's, 
this  variable  became  an  exogenous  policy  variable.   Endogenous  variables  in 
the  wholesale  market  include  production  of  gas  and  wholesale  prices  for  two 


sectors,  residential-commerciaL,  and  industrial.   These  are  derived  from 
a  price  markup  equation  linking  field  contract  prices  and  transportation 
costs  to  wholesale  prices  in  the  residential  and  industrial  sectors,  and 
from  demand  equations  in  the  two  sectors.   The  field  and  wholesale  markets 
operate  to  provide  reserves,  new  production,  and  wholesale  prices  for  new 
production,  once  ceiling  prices  for  reserves  are  set  by  F.P.C.  policy. 

The  model  is  estimated  by  pooling  cross-section  and  time  series  data 
from  the  1960 's  for  the  wholesale  market,  and  for  supply  in  the  field  market, 
over  regional  break-downs  that  are  particular  to  each  part  of  the  model. 
Our  method  for  pooling,  and  the  estimation  techniques  as  a  whole,  are 
treated  in  the  third  section.   Section  4  describes  the  statistical  results 
and  the  fifth  section  describes  preliminary  results  from  a  simulation  of 
the  model  for  different  price  increase  or  "deregulation"  policies  of  the 
Federal  Power  Commission. 

2.    Structure  of  the  Model 

The  model  is  a  policy  model,  since  relationships  describing  field  price 
regulation  are  introduced  and  serve  to  change  the  structure.   Both  the  well- 
head price  in  the  field  market  and     the  wholesale  price  of  gas  sold  by 
pipeline  companies  are  determined  in  a  different  way  by  regulation.  When 
regulation  does  not  interfere  with  the  equilibrium  of  supply  and  demand, 
prices  and  quantities  in  both  markets  are  endogenous  variables.   However, 
if  the  regulated  ceiling  price  is  lower  than  the  equilibrium  market  price, 
we  no  longer  observe  both  the  demand  curve  and  the  supply  curve  (depending 
on  the  market)  and  price  is  exogenous  in  one  or  both  sets  of  markets. 


In  the  field  market  the  supply  variable  is  the  change  in  reserves, 
AR,  that  can  be  "committed"  to  pipelines  and  local  industrial  buyers  for 
production  over  a  ten  to  twenty  year  period.   The  change  in  reserves  comes 
about  from  additions  in  both  non-associated  gas,  and  that  associated  with 
oil  reserves,  that  result  from  new  discoveries  and  extensions  and  revisions 
of  earlier  discoveries.    The  demand  variable  in  the  field  market  is  the 
wellhead  price  that  pipelines  are  willing  to  pay  for  new  reserves.   After 
1961,  however,  this  price  variable  became  the  exogenous  F.P.C.  ceiling 
price,  since  it  was  lower  than  the  price  that  would  have  resulted  from  a 
market  equilibrium  (cf . [4]  Chapter  4,  and  [10]   ). 

In  the  wholesale  market  the  demand  variable  is  production  of  additional 
gas  AQ  for  delivery  to  public  utilities  and  industrial  consumers  as  a 
function  of  wholesale  delivery  prices,  the  prices  of  alternative  fuels, 

and  economy-wide  variables  that  determine  the  size  of  markets  for  energy. 

2 
Because  the  wholesale  market  is  oligopolistic,   there  is  no  explicit  "supply 

variable,  but  rather  wholesale  prices  are  endogenous  and  determined  by 

pipeline  markups  on  the  field  price  and  the  marginal  costs  of  transmission. 


The  findings  are  "supply"  variables.   There  are  other  accounting  defi- 
nitions of  reserves  different  from  AR,  however,  that  are  important  for  as- 
sessing the  stock  available  at  any  one  time.   The  stock  is  equal  to  AR,  less 
subtractions  due  to  production  and  any  net  change  in  storage.   The  levels 
are  affected  by  the  existence  ijf  repressuring  (which  adds  to  reserves), 
field  use,  transportation  losses,  etc.   A  flow  diagram  for  non-associated 
reserve  levels  in  shown  in  Figure  1  (page  6) .   The  same  diagram  would  hold 
for  associated  reserves,  the  two  being  linked  at  Marketed  Production,  as 
shown  in  the  diagram. 

2 
There  are  actually  several  wholesale  "markets"  (just  as  there  are 

several  production  markets)  in  different  parts  of  the  country  that  are 

geographically  defined,  and  the  degree  of  monopoly  power  in  each  market 

is  different.   This  is  explained  in  more  detail  in  Section  3. 
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The  quantity  of  gas  sold  must  be  consistent  with  the  constraints 
imposed  by  existing  reserve  levels.   Pipeline  companies  can  only  deliver 
additional  gas  based  on  the  commitments  of  new  reserves  made  to  them  by 
producers.   These  commitments  are  usually  made  for  ten  to  twenty  year 
periods,  and  therefore  the  increment  of  gas  produced  to  meet  new  sales 
commitments  can  be  only  approximately  one-tenth  to  one-twentieth  of  the 
change  in  the  reser^/e  level.   If  demand  for  production  for  an  extended 
period  exceeds  this  allowable  level  of  extraction  from  new  reserves, 
rationing  has  to  take  place  (as  has  indeed  been  the  case  in  recent  years). 

A.   The  Field  Market 

Probably  the  sector  of  the  natural  gas  industry  most  difficult 
to  capture  in  a  conceptual  model  is  the  supply  of  new  reserves.   Most  of 
the  current  controversy  over  regulatory  policy  centers  on  this  sector  and 
whether  or  not  supply  has  been  too  low  as  a  result  of  past  regulatory  policy. 
Actual  additions  to  reserves  through  new  discoveries  are  realized  by  a 
complicated  process  involving  a  large  number  of  technological  factors,  and 
it  may  seem  naive  to  try  and  model  the  process  using  a  set  of  simple 
econometric  relationships.   Structural  equations  can  be  formulated,  however, 
that  do  link  several  economic  variables  that  are  important  in  gas  reserve 
additions  and  also  describe  most  simply  and  directly  the  regulatory  effects. 
Whether  or  not  the  econometric  description  is  adequate  can  be  judged  in 
part  from  the  logic  of  the  interrelationships  and  in  part  from  the  statis- 
tical results. 


A.l   Supply  of  Reserves 

The  major  component  of  additions  to  reserves  is  new  discoveries 
of  both  non-associated  and  associated  gas  (associated  includes  "dissolved" 
gas  recovered  from  oil  production,  as  well  as  "free"  natural  gas  forming  a 
cap  in  contact  with  crude  oil) . 

The  process  of  discovering  gas  begins  with,  and  is  driven  by,  the 
drilling  of  wells.   Discovery  is  a  process  for  obtaining  productive  inputr 
the  same  as  any  other  oil  and  gas  industry  process.   The  results  follow  when 
there  are  sales  of  new  reserves  AR;  to  obtain  these  reserves  at  least  costs, 
the  exploratory  firm  minimizes  costs  —  expenditures  on  labor  L  and  well 
drilling  W  —  subject  to  the  production  function  AR  =  f (L,W)  as  a  constraint. 
The  least-cost  combination  of  inputs  results  in  supply  AR  =  f (PG,PW,CW,X. ) 
where  PG  is  the  field  contract  price  for  gas,  PW  is  the  industry  wage  rate, 
CW  is  unit  costs  for  wells,  and  X   is  a  parameter  i  in  the  production 
function.    In  modeling  new  discoveries  specifically,  however,  we  draw  from 
cind  expand  on  the  work  of  Erickson  and  Spann  [6],  Khazzoom  [14]  and  MacAvoy 
[18]. 

There  is  a  distinction  between  the  supply  curve  of  "exploratory  effort" 
and  supply  curve  of  actual  new  discoveries  (cf.  Erickson  and  Spann).   Economic 
incentives  not  only  influence  the  amount  of  exploration  but  also  determine 


Assume,  for  simplicity,  that  exploratory  activity  can  be  described 
by  the  Cobb-Douglas  function  AR  =  aL^^W^,  so  that  the  minimum  cost  combina- 
tion of  L,W  is  found  by  the  first-order  conditions  for  minimizing  costs 
C  =  P  -L  +  C  *W  subject  to  AR  =  aL^'^.   The  level  of  minimum  marginal 

costs  from  these  conditions  is  3C/oAR  =  f(AR,  P  ,C  ,a,6,Y)  and  with  the 

supply  of  reserves  set  at  the  level  at  which  price  of  nev;  reserves  PG 
equals  marginal  costs  SC/3AR,  then  supply  AR  =  f(PC,Pj^,C  ,a,B,Y)'   Here 

X  =  a,  x^  =  3,  and  x_  =  y. 


Its  characteristics  since  an  increase  in  incentives  not  only  leads  to  more 
wildcat  drilling  but  also  drilling  of  poorer  quality  prospects.   Any  measure 
of  the  results  of  exploration  such  as  the  success  ratio,  SR  (ratio  of  pro- 
ductive to  total  wildcats)  or  the  average  size  of  discovery  S,  is  a  function 
of  both  the  distribution  of  prospects  found  in  nature  and  the  risk  pre- 
ferences of  the  "wildcatters."   Thus,  total  new  discoveries  depend  on  (a)  the 
number  of  wildcats  drilled,  (b)  the  success  ratio,  and  (c)  the  average  si2e 
of  discovery  per  successful  wildcat,  and  since  the  latter  two  factors  de- 
cline with  a  rise  in  economic  incentives,  the  response  of  new  discoveries 
to  such  inventives  will  be  less  than  the  response  of  wildcatting  itself. 
The  supply  curve  of  wildcat  drilling  can  be  formulated,  however,  in 
terms  close  to  the  general  characterization  of  the  exploratory  process. 
The  "production  function"  process  consists  of  using  wells  and  labor  to 
select  from  the  prospects  those  able  to  provide  lowest  cost  reserves  AR. 
The  prospects  set  the  parameters  of  the  production  function.   Since  data 
on  the  characteristics  of  prospects  is  not  available,  it  is  necessary  to 
use  surrogates  —  the  characteristics  of  drilling  in  the  preceding  year. 
The  valuesof  the  lagged  variables  W    ,  SR  _  ,  ^^_^  '  ^^  ^  region  are  indi- 
cative of  the  year's  prospects  in  that  region,  because  they  measure  the 
information  available  to  wildcatters  regarding  current  results  in  that  area. 
Drilling  depends  on  the  "incentive"  variables  —  costs  per  foot  drilled, 
DC,  as  well  as  the  price  of  oil,  PO,   and  the  price  of  gas  PG  they  expect 


As  Khazzoom  [14]  points  out,  the  price  of  oil  comes  in  through  two 
effects:   directionality,  i.e.,  gas  discovered  in  the  search  for  oil, 
as  well  as  discovery  of  associated  gas. 


10 

to  receive.   With  discovery  wells  designated  as  WXG  for  those  proving  out 

2 
new  fields  as  well  as  for  extensions  to  the  size  of  known  fields,   and  with 

the  exact  functional  forms  for  the  equation  yet  to  be  determined,  we  expect 

the  following: 

^^t.j  =  ^(^^t,j'^°t,j'°^t,j'^Vij'^t-i,j'^Vi,j>     •  (i> 

Here  t  is  the  time  index  and  j  is  the  district  index  (a  total  of  19  Texas 
Railroad  Commission  and  other  substate  or  state  districts  are  used  to  obtain  a 
cross  section,  as  discussed  in  Section  3) . 

The  supply  functions  for  reserves  are  formulated  to  show  the  results 
from  drilling.   In  the  case  of  non-associated  gas,  new  discoveries  (DN)  are 
a  function  of  the  number  of  gas  wells  drilled  for  new  discoveries,  but  also 
of  the  geological  productivity  of  drilling  in  that  region.   As  an  indicator 
of  these  geological  conditions,  we  have  formulated  the  variable  MM.  as  the 
average  level  of  discoveries  per  well  in  district  j  over  the  entire  time 
period  t.   When  this  level  has  been  large  over  the  period,  it  can  be  inferred 


These  in  turn  may  be  weighted  distributions  of  past  prices  they  did       J 
receive.  ' 

2 
In  general,  wells  are  classfied  according  to  intent  when  they  are 

drilled  or  by  the  result  after  drilling  is  completed  (or  abandoned) .   The 
classification  according  to  intent  by  the  American  Association  of  Petro- 
leum Geologists  (AAPG)  consists  of  drilling  for  new  fields  (new  field 
wildcats),  drilling  for  new  pools  (new  pool  wildcats,  shallower-pool  tests, 
deeper-pool  tests)  or  to  extend  presently  known  pools  (outposts) .   After 
drilling,  wells  are  classified  as  "dry  holes"  if  they  are  unsuccessful,  or 
as  new  field  or  new  pool-discovery  wells  if  a  new  field  or  pool  is  dis- 
covered or  extension  wells  if  they  extend  the  size  of  presently  known  pools. 
We  center  attention  on  new  field  extension  wells,  as  representative  of  those 
that  take  the  development  of  the  field  from  the  level  of  a  single  wildcat 
well  to  the  level  of  reasonably  known  or  "discovered"  reserves. 


11 

that  potential  reservoir  content  has  been  greater  there  than  elsewhere. 
New  discoveries  of  non-associated  gas,  then,  will  have  the  form: 

DNj.^^  =  f(WXG^^^,MM^,DN^_^^^)  .  (2) 

Depending  on  the  relative  prices  of  gas  and  oil,  associated  gas  can 
be  viewed  as  an  input  factor  for  producing  oil  or  as  a  separate  output. 
As  Khazzoom  indicates  [14],  if  the  price  of  oil  is  considerably  greater  than 
the  price  of  gas,  then  gas  is  to  some  extent  an  input  for  producing  oil 
(since  gas  pressure  forces  the  oil  to  the  surface,  and  the  higher  the  gas 
pressure,  the  higher  is  the  oil  extraction  rate).   If  oil  prices  were  to 
be  reduced  relative  to  gas  prices,  the  gas-to-oil  production  ratio  should 
increase.   In  many  regions,  however,  state  conservation  laws  constrain  the 
gas-to-oil  production  ratio  to  be  within  some  range  to  conserve  the  "gas- 
input"  resources  so  as  ultimately  to  produce  more  oil.   This  effectively  puts 
upper  bounds  on  production.   In  this  case,  the  supply  of  associated  gas  will 
be  linked  to  the  production,  QO,  of  oil.   On  the  whole,  we  would  expect  dis- 
coveries of  associated  gas  to  be  more  a  function  of  oil  production  but  to 
also  depend  somewhat  on  well  drilling  for  oil  that  produces  both  oil  and  gas: 

°\.j  =  ^('^Vi,j'^'^^t,j'^Vi,j)     •  (^) 

Extensions  and  revisions  for  both  non-associated  and  associated  gas  are 
determined  by  previous  discoveries,  previous  or  present  production  of  both 
gas  and  oil,  and  perhaps  even  prices  (if  prices  induce  further  activities 
beyond  drilling  that  add  to  reserves  in  these  catesories) .   Extensions  and 
revisions  for  associated  gas  specifically  should  depend  on  the  rate  of 


production  of  oil  (QO) ,  particularly  in  those  regions  that  constrain  the 
gas-to-oil  production  ratio.  Our  equations  for  extensions  and  revisions 
will  thus  be  of  the  form: 

(Non-associated)  ; 

^^t.j  =  ^(^'t.j'5^-:.J'°VlJ'"^^t-lJ^Vl.j^  (^> 

(Associated) 

^t.j  =  ^^^^t,j'Q°tj'j?^\,j'°Vij'^Vij'^^t-i,j^  ^'^ 

where  prices  are  indicators  of  previous  discoveries  (2AR   .  ,  ^^,._i  •  ^'^'^ 
WXG   ^  .)  and  previous  extensions  and  revisions  are  included.   For  accounting 
purposes,  the  net  change  in  reserves  is  identically  equal  to  total  new  dis- 
coveries plus  total  extensions  and  revisions  less  marketed  production  (Q) , 
losses  (L)  and  change  in  underground  storage  AUS: 

^t.j  =  °\j  ^  ^tj  +  %,j  +  ^t.j  -  QtJ  -  \,i    -   ^"^.j-   (6> 

The  cross-sectional  index  for  the  change  in  underground  storage,  AUS ,  cor- 
responds to  the  production  regions.   Pipeline  companies  often  store  gas 
during  the  off-peak  season  in  the  vicinity  of  major  consumers  for  eventual 
use  to  satisfy  peak  demand.   Gas  inventories  thus  serve  as  a  substitute  for 
peak-load  transmission  capacity.   Although  the  gas  is  stored  in  the  con- 
sumption areas  (and  the  amount  of  storage  is  cost-limited  by  the  extent  of 
geological  conditions  favoring  inground  insertion  of  gas  in  those  areas), 
an  increase  in  inventories  is  a  drain  on  the  reserves  of  the  production 
region  which  supply   the  gas.   We  model  inventory  change  (as  well  as  losses 


In  the  econometric  work  that  follows,  separate  eauations  are  fitted 
for  extensions  and  for  revisions.  Lz   should  be  noted  that  they  all  follow 
from  this  general  formulation. 
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L)  as  a  fixed  proportion  of  production  that  is  particular  to  the  production 
region.   For  a  monthly  or  quarterly  model  this  would  be  unrealistic  because 
of  the  seasonal  fluctuation  in  inventory  change  that  depends  largely  on  the 
geographical  and  climactic  conditions  of  the  consumption  area.   Since  our 
model  is  annual,  though,  these  seasonal  fluctuations  are  not  important  and 
it  is  only  the  long-term  trend  in  inventory  change  that  we  need  concern 
ourselves  with. 

A. 2  Demand  for  Reserves 

The  demand  for  new  reserves  by  pipelines  is  specified  by  an  equation 
for  the  wellhead  price,  PG.  Some  time  after  1961,  however,  the  regulated 
ceiling  price  (PC)  has  to  be  below  the  price  that  would  have  resulted  from 
a  supply-demand  equilibrium  since  otherwise  the  ceiling  would  not  prevail. 
When  ceilings  became  effective,  excess  demand  became  a  necessary  result  and 
we  no  longer  observe  the  demand  curve  after  that  point.   This  is  shown 
graphically  in  Figure  2a  (before  ceilings)  and  2b  (after  ceilings). 


Before  Ceilings 


After  Ceilings 


Figure  2a 


Figure  2b 
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Before  ceiling  price  controls,  the  endogenous  wellhead  price  in  region  j  is 
a  function  of  total  new  reserves  ZAR   .  (as  shown  in  Figure  2a) ,  the  change 
in  reserves  in  "j"  alone  (because  costs  of  transmission  from  j  vary  with  the 
amount  of  reserves  in  that  region) ,  the  average  mileage  (M)  between  the 
center  of  the  production  region  and  the  distribution  area  (because  costs  of 
transmission  vary  with  distance),  and  the  total  quantity  of  new  production 
from  the  region  AQ   ..   But  the  buying  markets  are  not  the  producing  regions 
"j."   The  districts  have  to  be  aggregated  to  regional  "market"  areas,  i,  each 
of  which  consists  of  a  grouping  of  the  production  or  supply  districts  j  which 
are  identified  by  the  state  conservation  commissions  and  the  F.P.C.  to  be 
geological  units.   The  wellhead  price,  indexed  as  PG   .,  in  each  of  the 

nineteen  production  districts  depends  upon  reserves  there  and  is  aggregated 

1 
over  i  purchase  markets,  as  given  by: 

PG^j  =  ^S=i^t,j'^t,j'"t,j'^^t,j'^^t-l,j^'  ""  ~   '^eili^g  y^^"^ 

(7) 

=  PC  ;   t  -  ceiling  year  . 

In  some  markets,  even  without  price  controls,  demands  would  not  be  ob- 
servable, however.   Just  after  World  War  II,  the  few  existing  pipeline  com- 
panies had  monopsony  power  in  gas  field  markets  (cf.  MacAvoy  [17],  as  summarized 
by  Erickson  and  Spann  [6]).   In  a  monopsony  market,  individual  buyers  have 
an  effect  upon  the  prices  at  which  they  buy,  and  as  a  result  equilibrium 
price  and  quantity  lie  on  the  supply  curve  but  not  at  the  intersection  with 
the  demand  curve  as  in  Figure  2a.   If  a  profit-maximizing  pipeline  is  a 
monopsony  buyer  of  gas,  it  will  equate  demand  with  its  marginal  expenditure 


This  is  to  say  that  prices  on  individual  wellhead  sales  at  j  are  a 
function  of  the  total  amount  sold  in  the  market  i  made  up  of  this  group  of 
regions  j . 


15 

on  gas  so  as  to  find  AR*  —  but  it  pays  no  more  than  the  supply  price  P* 
for  this  reduced  volume  (see  Figure  3  below).   Thus,  the  equilibrium  price 
and  quantity  lie  on  the  supply  curve  but  below  the  demand  cur'/e. 


Monopsony  Market 


ME  (marginal  expenditure  curve) 


After  the  early  1950 's,  however,  this  condition  changed  as  the  network  of 
pipelines  grew  extensively.   Pipeline  companies  no  longer  enjoyed  a  monop- 
sony market  in  gas  fields,  except  perhaps  in  supply  regions  for  gas  to  the 
West  Coast  in  which  El  Paso  Natural  Gas  Company  was  a  very  large  source  of 
demand.   But  to  the  extent  that  monopsony  prevails,  prices  are  determined 
exogenous  to  the  system.   Whether  control  is  administered  by  the  F.P.C.  or  by 
pipeline  companies,  prices  PC   .  =  PC   .  in  a  "ceiling  year." 


16 


B.   The  Wholesale  Market 

Pipeline  companies  that  contract  for  reserve  commitments  from  producers 
sell  gas  on  the  wholesale  market  directly  to  industry,  and  to  local  public 
utilities  for  resale  to  residential,  commercial  and  industrial  consumers. 
We  have  identified  wholesale  market  regions  across  the  continental  United 
States,  with  each  region  representing  a  geographical  area  in  which  from  one 
to  several  pipeline  companies  compete  for  consumer  sales.   This  regional 
break-down  is  different  from  those  of  the  field  discovery  and  field  purchase 
markets,  primarily  because  wholesale  buyers  have  access  to  pipelines 
traversing  large  geographical  areas  and  they  operate  in  "markets"  covering 
large  regions.   Pipelines  from  the  same  field  demand  markets  may  feed  into 
two  or  more  regional  wholesale  markets,  and  pipelines  from  several  field 
demand  markets  may  all  feed  into  a  single  regional  wholesale  market.  (This; 
regional  break-down  is  discussed  further  in  Section  3.) 

Demands  in  these  markets  can  be  classified.   There  are  three  broad 
uses  of  natural  gas  by  industry,  and  for  each  use  the  quality  required  of 
the  gas  (and  thus  the  price  the  buyers  are  willing  to  pay)  is  different. 
Gas  used  for  chemical  processes  should  be  of  extremely  pure  quality  and  may 
be  sufficiently  unique  to  that  process  that  there  are  few  substitutes.   A 
second  use  for  industrial  gas  is  for  boiler  fuel,  and  here  the  gas  need 
not  be  very  pure  and  competes  with  oil  and  coal.   The  third  (and  smallest) 
use  of  industrial  gas  is  for  electricity  generation,  and  here  too,  the 
quality  of  the  gas  need  not  be  very  high  nor  the  substitutability  low,  since 
this  is  again  for  boiler  use.   Contracts  for  industrial  gas  are  also  made 
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on  either  a  firm  or  an  interruptible  basis.   Firm  contracts  require  that  gas 
be  supplied  throughout  the  year  at  a  more  or  less  constant  flow  rate,  while 
interruptible  gas  is  supplied  only  in  the  off-peak  season  when  there  is  ex- 
cess capacity. 

At  the  first  stage  of  model  construction,  the  three  industrial  gas 
demands  are  not  classified  separately.   One  reason  for  consolidating  dif- 
ferent demands  is  that  it  is  difficult  to  obtain  data  on  industrial  gas 
sales  broken  down  by  class  of  user;  pipeline  companies  must  report  to  the 
F.P.C.  gas  sales  to  each  industrial  firm,  but  they  do  not  report  the  ulti- 
mate use  of  the  gas.   It  is  necessary  even  at  the  first  stage,  on  the  other 
hand,  to  be  able  to  separate  "interruptible"  from  "firm"  sales,  particularly 
since  the  proportions  of  interruptible  to  firm  sales  are  different  in  the 
different  regional  markets  (the  greater  part  of  the  firm  industrial  sales 
are  made  in  the  Southwest,  close  to  field  sources),  so  that  industrial  de- 
mand elasticities  are  not  homogeneous  across  different  regions.   For  now, 
this  classification  is  made  on  the  basis  of  treating  the  firm- interruptible 
ratio  as  one  source  of  cross-sectional  heterogeneity. 

During  most  of  the  period  before  1968,  it  would  appear  that  wholesale 
markets  cleared  and  there  was  no  excess  demand  for  gas  by  consumers.  The 
equilibrium  quantity  of  gas  sold  by  pipelines,  and  therefore  produced  from 
reserves,  is  determined  by  the  intersection  of  the  consumer  demand  curve 


After  about  1968,  wholesale  markets  no  longer  cleared  and  the  demand 
curve  is  no  longer  observable.   It  can,  however,  be  extrapolated  from  pre- 
1968  data  and  in  this  way  we  can  obtain  a  measurement  of  excess  demand. 
We  shall  return  to  this  point  later. 
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with  a  curve  that  essentially  specifies  the  pipelines'  mark-up  on  the  field 
price.   The  price  mark-up  is  a  function  of  field  purchase  and  regional 
transmission  costs,  and  of  the  "cocpetitiveness"  of  both  gas  sales  and  inter- 
fuel  substitution. 

B.l   Demand  for  Gas 

In  this  model  the  demand  for  wholesale  gas  will  be  represented  by 
equations  that  relate  the  quantity  demanded  to  wholesale  price,  the  price 
of  alternative  fuels,  and  "market  size"  variables  such  as  population,  income, 
and  investment  that  determine  the  number  of  potential  consumers.   Two  demand 
equations  are  estimated,  one  for  pipeline  sales  to  retail  utilities  and  the 
other  for  pipelines  sales  to  direct  industrial  buyers.   In  both  cases,  the 
dependent  variable  is  new  demand,  6Q,  rather  than  the  total  level  of  demand. 
In  the  short  run,  as  Balestra  has  shown  for  residential  gas  [2],  the  level 
of  demand  should  be  relatively  price  inelastic  and  would  depend  on  stock 
variables  that  do  not  change  much  in  time  (e.g. ,  the  total  stock  of  gas- 
burning  appliances  for  residential  gas).   New  demand,  however,  should  re- 
spond to  the  price  of  gas  and  to  the  price  of  competing  fuels  (decisions  to 
buy  new  appliances,  for  example,  are  affected  by  fuel  prices).   The  new 
demand  for  gas,  6Q,  is  made  up  of  the  increment  in  gas  consumption  AQ  =  Q  -Q   •, 
and  of  replacement  for  continuation  of  old  consumption.   In  any  period, 
total  residential  and  commercial  gas  demand  could  be  considered  to  be  a 
function  of  the  stock  of  gas-burning  appliances.  A: 

Q^  =  A  •  A^  (8) 

where  X   is  the  (constant)  utilization  rate.   Then,  if  r  is  the  average  rate 
at  which  the  stock,  of  appliances  depreciates,  the  new  demand  for  gas  includes 
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rAA  ^^  .   The  total  new  demand  is   AQ  +  rXA  ^^ ,  or 

6Q^  =  AAj.  -  (l-r)XA^  .^  .  (9) 

Now  substituting    (8)    into    (9)    gives   us: 

6Qj.  =  Qj.   -    (l-r)Q^_^  (10) 

or 

<5Qt  =  AQj.  +  rQ^_^      .  (10a) 

Thus  the  new  demand  for  gas  is  the  sum  of  the  incremental  change  in  total 
gas  consumption  (AQ  )  plus  the  demand  resulting  from  the  replacement  of  old 
appliances. 

Our  _a  priori  assumption  is  that  new  demand  depends  on  prices  and  total 
income  (through  purchases  of  new  appliances),  but  the  level  of  total  demand 
is  itself  a  function  of  income  and  population.   Thus,  we  have  for  residential 
and  commercial  demand: 

^QS\,k  =  ^(^S\,k'P^,k'^,k'^\,k'^\,k'^QS\-l,k)  ^1^> 

where  PSR  is  the  wholesale  price  of  residential  and  commercial  gas,  PF  is 
a  price  index  of  competing  fuels,  Y  is  disposable  income,  N  is  population, 
in  regions  k,  and 

^^,k  =  ^^,k^^Vi,k  (^2) 

and 

<^\.k  =  ^t,k  ^  ^Vi,k   •  (^3> 

Balestra  [2]  distinguishes  between  two  depreciation  rates,  one  for  gas 
appliances  and  the  other  for  alternative  fuel-burning  appliances,  since 
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lifetime  for  appliances  using  alternative  fuels  may  be  different.   He  then 
estimates  the  two  depreciation  rates  by  estimating  an  equation  of  the  form: 

QSR^  =  a^  +  a^PSR^  -f  a^AN^  +  a3N^_^  +  a^AY^_^  +  a3Y^_^ 


(lA) 


cannot  be  justified  theoretically.)   The  all-fuel  depreciation  rate  comes 
out  of  equation  (14)  as  either  the  ratio  a./a„  or  a^/a, .   Thus,  the  equation 
is  over-identified,  and  the  depreciation  rate  can  be  obtained  only  by  esti- 
mating (14)  subject  to  the  constraint  of  a„/a„  =  a.   / a...       (The  resulting 
estimation  problem  is  non-linear,  but  Balestra  uses  an  iterative  method 
suggested  by  Houthakker  and  Taylor  [13]  to  obtain  an  estimated  depreciation 
rate  equal  to  0.11.)   The  two  depreciation  rates  needed  here  will  be  esti- 
mates in  keeping  with  these  procedures. 

The  demands  of  direct  industrial  purchasers  should  not  differ  greatly 
from  those  of  utility  purchasers  shown  in  (11).   But  there  are  some  distinc- 
tive variables.   As  mentioned  above,  sales  of  industrial  gas  are  made  on 
either  a  firm  or  interruptible  basis,  with  the  price  for  interruptible 
contracts  considerably  lower  than  that  for  firm  contracts.   Because  in  this 
model  we  do  not  disaggregate  firm  from  interruptible  sales  and  because  the 
ratio  of  firm  to  interruptible  sales  is  different  in  different  parts  of  the 
country,  we  will  introduce  into  our  demand  equation  for  industrial  gas 
an  exogenous  variable,  FIR    ,  for  the  ratio  of  firm  to  interruptible 
sales  of  industrial  gas  during  year  t  in  wholesale  market  region  k.   Since 
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firm  and  interruptible  prices,  we  would  expect  the  variable  FIR    to 
appear  in  the  demand  equation  with  a  positive  coefficient. 

Although  some  industrial  uses  of  gas  require  a  high  level  of  chemical 
purity,  on  the  average  the  cleanliness  and  ease  of  handling  of  industrial 
gas  are  not  as  important  as  in  residential  and  commercial  consumption.   As 
a  result,  industrial  gas  can  usually  be  sold  only  at  prices  equivalent  to 
or  lower  than  the  price  of  competing  fuels.   We  would  expect,  then,  that 
the  exogenous  price  index  of  competing  fuels,  PF,  will  be  an  important 
variable  in  the  industrial  wholesale  demand  equation. 

As  was  the  case  for  the  residential  and  commercial  sector,  our  demand 
equation  for  industrial  gas  will  use  new  demand,  <5QI,  as  the  dependent 
variable.   We  cannot  expect  the  depreciation  rate,  r  ,  however,  to  be  the 
same  as  it  was  in  the  residential  and  commercial  sector.   Total  industrial 
demand  depends  also  on  market  size  variables,  such  as  population  and  income. 
But,  since  new  gas  in  the  industrial  sector  is  for  use  in  producing  more 
output,  it  might  be  expected  that  industrial  capital  investment  would  be 
the  most  important  determinant  of  new  energy  demands.   Our  demand  equation 
will  be  of  the  form: 


6QI   ,  =  f (PI   ,  ,PF   ,  ,FIR   ,  ,6^N   ,  ,6^Y   ,  ,K   ,  ,6QI   ,  , 
^  t,k       t,k'   t,k'    t,k'  I  t,k'  I  t,k'  t,k'  ^  t-l,k 


)       (15) 


with      «QIt,k  =  ^Q^t,k+^lQVl,k  (16) 

^I^.k  =  ^^t.k  ^  ^lQ^t-l,k  <^7> 

Vt.k^n.k  +  '^iVi.k  (^«> 

and  K    is  corporate  capital  investment  in  wholesale  market  k  in  year  t. 
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B.2  Gas  PrlcinR  at  Wholesale 

With  large  numbers  of  connecting  pipelines  between  producing  dis- 
tricts and  consumers  developed  during  the  1960's,  the  industrial  and  utility 
buyers  in  the  wholesale  markets  were  faced  with  a  variety  of  "supply"  con- 
ditions.  In  some  markets,  there  were  one  or  two  sources  of  new  gas,  so  that 
prices  were  fixed  by  the  pipelines  to  meet  demands.   That  is,  strictly 
speaking,  there  was  no  "supply  function,"  but  rather  a  price  mark-up  on 
field  purchases  set  by  costs  and  Federal  Power  Commission  regulation.   In 
other  markets,  there  were  more  transporters  of  gas  and  the  mark-ups  were 
as  limited  by  the  conditions  of  inter-pipeline  rivalry  as  by  those  condi- 
tions set  by  regulatory  procedures.   But  even  under  extensive  inter-pipeline 
competition,  the  process  is  the  same:   the  wholesale  prices  are  set  by  th* 
pipelines  equal  to  the  wellhead  price  plus  a  mark-up  proportional  to  the 
marginal  costs  of  transmission  and  distribution.   The  "proportionality" 
factor  varies  by  the  extent  of  the  rivalry,  and  the  nature  of  regulation. 

The  marginal  costs  of  transmission  can  be  found,  at  least  through  an 
"averaging  process"  similar  to  that  in  computing  regression  equations. 
We  assume  to  begin  with  that  in  an  existing  pipeline  system  average 
variable  costs  of  transmission  are  constant  up  to  some  capacity,  and  then 
rise  with  infinite  slope.   This  lets  us  net  out  the  marginal  cost  of  pur- 
chasing gas  at  the  wellhead,  as  long  as  pipelines'  contracts  for  new 
reserves  can  satisfy  final  demand  for  new  production.   Thus,  there  are 
two  factors  that  can  result  in  infinite  marginal  costs:   limited  trans- 
mission capacity,  and  a  limitation  on  the  ability  of  pipeline   companies 
to  deliver   enough  gas  production  from  the  available  new  contracts  for 


23 

field  reserves.   We  take  the  marginal  costs  of  transmission  to  be  a 

function  of  mileage  M,  ,  capacity  of  the  lines  transmitting  from  region 

k  that  year  v  ,  ,  and  the  ratio  of  firm  to  interruptible  sales  FIR  ,  . 
t,k  t,k 

Initially,  we  will  estimate  some  very  simple  wholesale  price  equations 
of  the  form: 


and 


^^^t.k  -  ^^t,k  =  ^^.k'\,k'^^\.k)  (1^> 


^^t.k  -^^t.k=  ^^.k^^^.k'^'^.k)   •  ^'"^ 


wellhead  prices  from  those  production  regions  which  feed  pipelines  that 
transmit  gas  to  that  wholesale  market  k.   The  mileage  variable  tl  and 
capacity  variable  v  ,  are  similarly  defined,  i.e.,  in  terms  of  the  average 
values  of  these  variables  from  production  regions  to  the  wholesale  market 
region  k.   The  "mark-up"  is  thus  the  excess  of  the  price  over  field  costs 
allowed  by  competition  and  the  Commission,  but  as  shown  by  the  size  of 
coefficients  of  the  independent  variables  in  marginal  transmission  costs. 
The  "determined"  price  is  approximated  by  the  regression  equation  coeffi- 
cients. 

When  the  wholesale  price  is  regulated  to  be  a  mark-up  over  the  well- 
head price  and  the  marginal  cost  of  transmission  (as  it  is  in  the  case  of 
residential  and  commercial  gas),  excess  demand  can  exist 


This  is  in  the  tradition  of  models  of  pipeline  costs  such  as  that  of 
Wellisz  [20],  V7here  transmission  capacity  is  tha  limiting  factor.   Further- 
more, in  recent  years,  when  there  was  excess  demand  in  wholesale  markets, 
it  was  new  reserves  that  placed  a  lim.itation  on  wholesale  production  capacity. 
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in  the  wholesale  market.   This  is  ihown  in  Figure  4.   When  demand  increases 
from  DD  to  D'd',  the  wholesale  price  is  regulated  to  be , a  mark-up  over 
marginal  costs  (MC^.,  the  marginal  cost  of  purchase  in  the  field,  is  just 
the  wellhead  price,  PG) ,  and  marginal  costs  become  infinite  at  AQ    ,  then 
there  will  be  excess  demand  as  shown.   The  question  is,  at  what  point 
AQ    does  the  marginal  cost  curve  become  vertical?   This  occurs  at  the 
point  at  which  new  production  is  constrained  by  the  lack  of  productivity 
from  the  given  level  of  new  reserves.   Generally,  we  would  expect  to  have 


Wholesale  Price 


Field  price 


MC  +  MC^  =  MC  +  PG 


Figure  4 
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OAR  =  AQ^^  (21) 

with  a  set  at  a  limit  by  the  full  flow  characteristics  of  the  new  reservoir 
put  under  contract.   The  problem  is  that  it  is  difficult  to  use  past  data 
to  measure  a  since  up  until  about  1968,  at  least,  production  demands  were 
such  that  AQ  <  AO    .   We  will  attempt,  however,  to  use  cross-sectional 
data  over  production  regions  in  1969  and  1970  to  run  regressions  of  the 
form 

AQj.^j  =  a  AR^^j  +  BRj.^j  (22) 

under  the  assumption  that  in  those  years  AQ  =  AQ    (we  will  exclude  those 

districts  for  which  this  was  not  believed  to  be  the  case) . 

Once  we  have  a  measure  of  AQ_   ,  we  can  extraoolate  the  demand  curves 
Tnax 

estimated  using  pre-1968  data  (since  equilibrium  production  and  price  did 
not  lie  on  these  curves  after  1968)  to  determine  excess  demand  in  recent 
years.   In  other  words,  we  can  simulate  the  model  past  1968  (thus  the 
demand  curve  D'd'  is  just  the  demand  curve  DD  simulated  into  the  future) 
using  a  constraint  on  production: 

Z  AQSR^  ^  +  L   AQI^  ^   =  AQ^^""   ;  t  -  1968     .  (23) 

k      '     k     ' 

This,  unfortunately,  does  not  take  into  account  the  complicated 

interconnections  between  regional  field  markets  and  regional  wholesale 

markets.   Our  constraint 

E  AQ^  a  Z  AR  (24) 

k   ^     j   J 

is  actually  somewhat   incorrect.   In  fact,  there  should  be  a  separate  con- 
straint for  each  wholesale  market  region  k: 
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l,k 
where  the  above  summation  is  taken  over  those  i's  that  feed  into  the 
particular  k,  and  3.  ,  is  the  fraction  of  AR.  going  to  region  k. 


3.    The  Use  of  Cross-Section  and  Time-Series  Data  for  Estimation 

The  equations  of  this  model  describe  the  dynamic  equilibria  of  two 
distinct  (but  interrelated)  groups  of  regional  markets,  namely,  field 
markets  and  wholesale  markets.   The  supply  side  of  each  regional  field 
market  is  determined  by  the  discovery  and  accumulation  of  new  reserves 
in  several  different  production  districts.   In  fact,  approximately  20 
production  districts  determine  the  supply  of  gas  reserves  in  6  or  so 
regional  field  markets.   The  demands  for  gas  on  the  wholesale  level  are 
determined  in  about  5  different  regional  markets.   As  well,  the  field 
and  wholesale  markets  are  themselves  interrelated.   Thus  the  cross-sectional 
break-down  will  be  different  in  the  different  sectors  of  the  model. 

The  time-bounds  over  which  equations  can  be  estimated  will  also  be 
different  for  different  equations.   This  is  the  case  for  two  reasons. 
First,  the  structure  of  the  models  changes  over  time.   To  begin  with,  the 
degree  of  monopsony  and  monopoly  power  held  by  pipelines  in  different  re- 
gional markets  has  been  changing  in  time.   But  even  more  fundamental  is 
the  fact  that  when  regulation  becomes  effective  (i.e.,  when  the  ceiling' 
price  is  below  the  equilibrium  price)  the  demand  curve  is  no  longer  ob- 
servable and  new  reserves  are  defined  by  the  supply  curve  and  the  ceiling 
price.   The  second  reason  for  the  use  of  different  time-bounds  is  that  the 
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availability  of  statistical  data  is  different  for  different  variables. 
For  many  disaggregated  variables,  particularly  related  to  field  supply, 
data  has  been  collected  and  compiled  only  after  1965;  this  imposes  con- 
straints on  the  tine-bounds  over  which  we  can  estimate  particular  sectors 
of  the  model,  but  not  other  sectors.   The  equations  of  this  model,  then, 
will  be  estimated  by  pooling  cross-section  and  time-series  data,  but  usinf 
cross-sections  and  time-bounds  that  are  different  for  different  parts  of 
the  model. 

The  differences  in  coverage  and  time-lengths  is  extensive.   Table  1 
contains  a  sunmiary  of  the  model's  equations  (identified  by  dependent 
variable)  grouped  according  to  cross-sectional  break-down  and  time-bounds. 
The  remainder  of  this  section  will  explain  the  geographical  break-down  in 
more  detail,  and  will  also  discuss  some  of  the  statistical  aspects  of 
pooling  heterogeneous  cross-sections  and  time-series. 

A.    Identification  of  Regional  Markets 

Although  data  are  available  for  new  reserves  and  field  price  for 
production  districts,  these  districts  do  not  by  themselves  constitute  a 
"market"  in  the  usual  sense  of  the  word,  since  many  districts  may  be  selling 
gas  to  the  same  group  of  pipeline  companies.   We  define  a  "field  market" 
as  a  distinct  geographical  region  i  where  gas  producers  in  all  of  the  pro- 
ducing regions  j  in  competition  with  each  other  sell  reserve  commitments 
to  pipeline  companies  also  in  competition  with  each  other  as  buyers. 

As  one  would  expect,  the  identification  of  regional  field  markets  is 
not  a  straightforward  matter.   The  AGA  Reserve  Committee  has  identified  ten 
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Table  1: 

THE  MODEL'S  STOCHASTIC  EQUATIONS  GROUPED  BY 
CROSS-SECTIONAL  BREAK-DOWN  AND  TIME-BOUNDS 


1.  Cross-Section  Break-do^m  by  Production  District  (j's) 

Equations;  Time  Bounds 

Wells  (WXG)  1954  -  1971 

Additions  to  Reserves  (DN,  DA,  XN,  XA,  RN  1966  -  1971 
RA) 

2.  Cross-Section  Break-down  by  Field  Markets  (I's) 


Wellhead  Price  (PG)  1954  to  date  ceiling 

price  is  effective 
(approximately  1961) 


3.    Cross-Section  Break-down  by  Wholesale  Market  (k's) 

Wholesale  Demand  (QRC,  QI)  1954  -  1971 

Wholesale  Price  (PRC,  PI)  1954  -  1971 


29 

regional  field  markets,  but  these  markets  do  not  seem  to  be  distinct, 
i.e.,  there  are  reserve  commitments  across  markets.   Therefore,  we  make 
an  identification  of  six  field  markets  based  on  geography,  but  also  on 
judgments  about  where  competition  between  districts  does  and  does  not 
exist.   The  two  alternative  classifications  of  field  markets  are  sum- 
marized in  Table  2.   Estimation  of  the  model's  field  market  demand  equa- 
tions can  be  made  using  both  classifications,  although  on  the  first  round 
of  construction  only  the  M-P  classifications  are  used;  eventually,  con- 
trasting  results  will  be  used  to  help  choose  the  most  reasonable  classi- 
fication. 

Regional  wholesale  markets  are  defined  and  identified  in  much  the 
same  way.   We  define  a  "wholesale  market"  as  a  distinct  geographical  region 
where  pipeline  companies  in  competition  with  each  other  sell  gas  to  public 
utilities  and  industrial  consumers  who  also  are  in  competition  with  each 
other  as  buyers.   The  identification  of  these  markets,  also  not  straightfor- 
ward, is   based  -on  an  examination  of  pipeline  locations  and  regional  data 
on  sales  of  gas.   Our  break-down  of  wholesale  markets  is  shown  in  Table  3. 


30 

Table  2; 
ALTERNATIVE  CLASSIFICATIONS  OF  FIELD  MARKETS 


AGA 


MacAVOY/PINDYCK 


1.  Appalachian 

Kentucky 

Maryland 

New  York 

Ohio,  Pennsylvanl£ 

Tennessee 

Virginia 

West  Virginia 

2.  Southeastern 

Alabama 
Florida 
Mississippi 

3.  South  Central  Area 

Arkansas 
Louisiana  North 


Field  Market  2 

Louisiana  Southeast 

Texas  1 

Texas  2 

Texas  3 

Texas  4 

Texas  5  &  6 


South  Louisiana  Area 

Louisiana  South 
Louisiana  Offshore 


Field  Market  3 

Louisiana  South 
Louisiana  Offshore 


5.   Texas  Gulf  Coast  Area 

Texas  1 
Texas  2 
Texas  3 
Texas  4 


This  classification  excludes  states  without  significant 
additional  reserves  in  the  1960's. 


31 


Table  2  (continued) 


Northeast  Texas  Area 

Texas  5 
Texas  6 
Texas  7B 
Texas  9 


MacAVOY/PINDYCK 


West  Texas-Southeast 
New  Mexico  Area 

Texas  7C 

Texas  8 

Texas  8A 

New  Mexico  Southeast 


Field  Market  4 


Texas  7C  &  8 

New  Mexico  Southeast 


Mid-Continent  Area 

Illinois 

Indiana 

Iowa 

Kansas 

Oklahoma 

Michigan 

Missouri 

Texas  10 


Field  Market  1 

Kansas 

Oklahoma 

Texas  10 

Texas  7B,  8A  &  9 


9.    Rocky  Mountain  Area 

Colorado 
Montana 
Nebraska 
.  New  Mexico  Northwest 
North  Dakota 


Field  Market  6 

Colorado  &  Utah 

Montana 

New  Mexico  Northwest 

Wyoming 


10.   Pacific  Coast  Area 

Alaska 

Arizona 

California 

Oregon 

Washington 


Field  Market  5 
California 
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Table  3:   District  Breakdown  for  Wholesale  Markets 


A.      SALES  FOR  RESALE 


1. 

Alabama 

2. 

Arizona 

3. 

Arkansas,  Missouri 

4. 

California 

5. 

Colorado,  New  Mexico,  Utah 

(Rocky  Mountain) 

6. 

Connecticut,  Massachusetts 

Rhode  Island  (New  England) 

7. 

Florida 

8. 

Georgie 

9. 

Idaho,  Wyoming 

10. 

Indiana 

11. 

Illinois 

12. 

Iowa 

13. 

Kansas 

14. 

Kentucky 

15. 

Louisiana 

16. 

Maryland,  Virginia 

17. 

Michigan 

18. 

Minnesota 

19. 

Mississippi 

20. 

Nebraska 

21. 

Nevada,  Oregon,  Washington 

22. 

New  Jersey 

23. 

New  York 

24. 

Ohio 

25. 

Oklahoma 

26. 

Pennsylvania 

27. 

Tennessee 

28. 

Texas 

29. 

West  Virginia 

30. 

Wisconsin 

B.   MAINLINE  SALES 


1. 

Alabama,  Georgia 

2. 

Arizona 

3. 

Arkansas 

4. 

Colorado,  Wyoming 

5. 

Florida 

6. 

Illinois 

7. 

Iowa,  Nebraska 

8. 

Kansas 

9. 


Kentucky,  Ohio,  Pennsylvania, 
West  Virginia  (Appalachian) 


10. 

Louisiana 

11. 

Minnesota 

12. 

Mississippi 

13. 

Missouri 

14. 

Oklahoma,  Texas 

15. 

Tennessee 
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B.    Estimation  Using  Pooled  Data 

The  cross-sectional  breakdowns  described  above  were  chosen  in  order 
to  identify  functional  economic  markets,  and  no  attempt  was  made  to  insure 
that  the  resulting  districts  will  be  homogeneous.   It  is  very  likely,  in 
fact,  that  observations  across  districts  are  not  homogeneous.   The  whole- 
sale market  in  the  northeastern  U.S.  probably  has  different  characteristics 
from  that  in  the  southeast  due  to  different  weather  conditions  and  different 
degrees  of  industrialization,  and  gas  discovery  rates  in  one  Texas  Railroad 
Commission  district  may  be  different  from  those  in  another  due  to  geological 
differences.   If  district  heterogeneity  is  not  accounted  for  when  specifying 
the  equations  of  our  model,  we  can  expect  the  additive  error  terms  to  be 
autocorrelated  across  districts,  just  as  time  trends  that  are  not  accounted 
for  will  result  in  error  terms  that  are  autocorrelated  across  time. 

This  section  will  briefly  review  some  of  the  standard  methods  in 
dealing  with  heterogeneity  when  estimating  with  pooled  data  and  will  then 
suggest  an  approach  to  be  used  in  this  study. 

B.l  Problems  of  Autocorrelated  Errors 
We  are  estimating  equations  of  the  form 

Y.^  =  Xj^3  +  e.^  (26) 

where  j  is  a  cross-section  index  (j=l,...,N)  and  t  is  the  time  index  (t=l,...,T) 
If  the  equation  is  estimated  by  ordinary  least  squares  (OLS),  the  Durbin- 
Watson  statistic  might  indicate  the  presence  of  autocorrelation  in  the  error 
terms,  but  it  will  not  tell  us  what  part  of  the  autocorrelation  is  across 
time  and  what  part  is  between  cross-sections.   Furthermore,  the  standard 
correction  techniques,  such  as  Hildreth-Lu  [11],  cannot  be  used  since  the 
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autocorrelation  is  two-dimensional.   If  the  errors  are  autocorrelated  and 
OLS  is  used,  we  can  expect  that  the  resulting  estimates  will  at  best  be 
consistent  and  unbiased,  but  inefficient.   If,  however,  the  equation  also 
contains  lagged  dependent  variabJ es  (e.g.,  if  Y    -  is  one  of  the  components 
of  X   ),  or  independent  variables  referenced  across  districts,  then  the 
estimates  will  also  be  inconsistent  [lO]. 

The  problem  of  autocorrelation  in  the  cross-section  dimension  is  often 
the  result  of  a  mis-specification  that  can  be  anticipated.   Suppose,  for 
example,  that  new  discoveries  of  gas  (ND)  is  believed  to  be  linearly  re- 
lated to  the  number  of  wells  drilled  (W) ,  so  that  the  equation  to  be  esti- 
mated is 

^j,t=^0-^  Vj,t-^^j,t    •  (27) 

It  is  reasonable,  however,  to  believe  that  geological  differences  will  make 

some  regions  richer  in  gas  than  others,  and  therefore  the  wells  in  those 

regions  will  have  a  higher  average  "output."  Perhaps  in  any  given  year, 

the  same  number  of  wells  per  district  in  each  of  two  different  districts 

j  and  j'  can  be  expected  to  result  in  different  amounts  of  discoveries. 

It  is  easy  to  see  that  this  situation  will  result  in  autocorrelated  errors 

in  equation  (27) . 

Consider  two  different  districts,  j  and  j',  with  average  "output  ratios" 

given  by 

T  ND 

Y  ^  (-1^)   =  a  (28) 

t=l  "j,t     J 

and  ^  I    ("IT"^)  =  «i'         •  (29) 

^  t=l  "j',t     J 
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Thus,  if  the  number  of  wells  in  these  two  districts  were  always  the  same, 
we  would  still  expect  to  find  on  the  average  that 

A  model,  then,  that  would  account  only  for  the  geological  differences  between 
districts  j  and  j'  would  be 

where  the  error  term  e*        is  independent  of  j .   Now  if  this  equation  (31)  is 
substituted  for  ND .    in  (27),  and  the  resulting  equation  is  written  with 
e    on  the  left-hand  side,  we  have 

But  ND.,    =  Bp.  +  6,W.,    +  e.,  ^,  and  substituting  this  into  (32)  gives 
J>t    (J    Ijjt     J>t 


^i,t  -  '3j'''i"r.t  -  Vj..  +  ^jj'^o  -  ^0  -^  =j,t  -^  Oiy'y.t    <"' 


^^'i.^'y.^^-'iA  ""> 


Errors  autocorrelated  in  time  can  occur  in  the  same  way.   Consider 
the  regression  equation  Y  =  3X  +  e  with  an  unexplained  time  trend; 

e.g.,  Y^  =  pY|._^  and  X^.  =  pX^_^.   Then,  e^_^  =  Y^_^  -  BXj._^  =  pY^  -  6pX^  =  pe^., 

so  that  E[e^  e^_^]   =  pa^. 
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B.2  Methods  of  Estimating  with  Pooled  Data 

Several  methods  exist  for  dealing  with  cross-sectional  heterogeneity 
when  estimating  equation?  using  pooled  data.   Probably  the  two  most  common 
approaches  are  to  introduce  district  dummy  variables,  or  else  to  use  a 
lagged  dependent  variable  to  account  for  district  heterogeneity. 

District  dummy  variables  have  been  used  in  natural  gas  studies  both 
by  Balestra  [  2]  and  by  Khazzoom  [14]  although  in  both  cases  the  statistical 
significance  of  the  dummy  variables  was  sometimes  questionable.   There  are 
two  basic  problems  with  district  dummies,  however,  that  would  make  their 
use  in  this  study  unappealing.   The  first  is  that  they  use  up  a  large  number 
of  degrees  of  freedom.   Some  of  the  equations  of  this  model  are  estimated  by 
pooling  data  over  some  20  districts,  but  over  only  5  or  6  years,  and  the 
introduction  of  district  dummies  would  thus  use  some  20%  of  the  degrees 
of  freedom.   The  second  problem  is  that  district  dummies  do  not  really  explain 
in  any  kind  of  systematic  way.   As  we  will  see,  heterogeneity  in  an  equation 
such  as  (27)  can  better  be  dealt  with  by  introducing  a  constructed  variable 
that  explains  past  differences  between  districts.  ' 

The  use  of  the  lagged  dependent  variable  can  be  a  convenient  way  to 
deal  with  heterogeneity,  particularly  if  the  errors  are  not  autocorrelated 
in  time.   If  the  equation  already  contains  lagged  independent  variables, 
then  the  use  of  the  lagged  dependent  variable  will  not  use  up  any  extra 
degrees  of  freedom.   One  problem  that  can  arise,  however,  is  an  estimated 
value  greater  than  one  for  the  coefficient  of  the  lagged  dependent  variable. 
Aside  from  the  difficulty  of  interpreting  such  a  result,  the  estimated 
equation  will  be  dynamically  unstable.   Balestra  had  this  problem  in  his 
estimates  of  demand  for  residential  gas.   He  frequently  obtained  values  for 
the  coefficient  of  the  lagged  dependent  term  that  were  greater  than  one, 
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implying  a  negative  depreciation  rate  for  gas-burning  appliances. 

Another  approach  is  to  assume  that  the  error  terms  are  made  up  of 
components  that  originate  from  different  sources  and  that  therefore  have 
different  variances,  and  then  gain  efficiency  in  the  estimates  by  using 
a  variation  of  generalized  least  squares  (GLS).   The  "residual"  or  "error 
components"  model  was  first  suggested  by  Kuh  [15],  and  later  generalized 
and  applied  by  Balestra  and  Nerlove  [  3  ]  and  Wallace  and  Hussain  [ 19] . 
Basically,  this  model  assumes  that  the  error  term  of  equation  (27)  is 
made  up  of  three  independent  components,  one  of  which  is  associated  with 
time,  one  with  the  cross-sections,  and  the  last  an  independent  random 
variable  across  both  time  and  cross-sections,  i.e.,  e.    is  given  by 


^j.t  =  ^j  -'^■'"jt 


(35) 


2   2       2 
with  variances  o    ,   o    ,    and  a  .   It  is  assumed  that  u.,  v^,  and  w_  are  all 
u    V        w^  J    t        J  t 

independent  of  each  other  and  that 

E[u.u.,]  =  0   for  j  ?*  j '  ,  E[v  v  ,  ]  =  0  for  t  ?*  t',  and 

E[Wjj.Wj,^]  =  E[w^^,Wj^]  =  E[w  ,^,w^^]  =  0  for  j  ji  j'  and  t  ?^  t'. 

Given  these  assumptions  about  the  error  vector  £,  ,  we  can  write  its 
variance-covariance  matrix  as 


R  =  E[ee']   =  o'^A  +  o^B  +  a^l„^ 
u     V     w  NT 


(36) 


Note  that  Q.   is  an  NT^NT  matrix, 
and  B  are  NTXNT  matrices  defined  by 


T 
0   J 


0  .  .  .  0 

„  .  .  .  0 


(37) 
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.   In 


(38) 


.    In 


2   2   2 
If  the  variance  components  a    ,    o    ,   o     are  known,  then  the  minimum  variance 


estimate  of  6  is  given  by  the  GLS  estimate 
-1   -1   -1 

e  =  (x'r2  -^x)    x'f2  -"y  . 


(39) 


If  the  variance  components  are  not  known  (which  would  presumably  be  the 
case),  then  Zellner's  method  [21]  can  be  used,  where  consistent  (but  ineffi- 
cient) estimates  of  6  are  obtained  by  OLS,  the  residuals  are  used  to  obtain 

2    2       2 
consistent  estimates  of  a  ,0  ,  and  o    ,  and  GLS  is  finally  used  to  obtain  a 
w    u       V 

new  (and  efficient)  estimate  of  g. 

The  problem  with  this  method  is  that  while  it  accounts  for  differences 
in  the  variances  of  the  error  components,  it  does  not  account  for  hetero- 
scedasticity  or  autocorrelations  within  each  error  component.   Thus,  if  the 
error  component  that  is  cross-sectionally  generated  is  itself  heteroscedastic 
or  if  its  elements  are  autocorrelated,  we  will  still  obtain  inefficient 
estimates  for  B  (although  the  estimates  will  be  more  efficient  than  those 
generated  by  OLS).   Also,  the  method  is  computationally  expensive,  since 

the  variance-covariance  matrix  fi  is  a  function  of  two  parameters  which  must 

2   2  2   2 

be  "scanned,"  namely  the  variance  ratios  p,  =  O   /o      and  p„  =  0  /a  . 

'       ■'  1    u  w      2    v  w 
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Estimation  of  the  equations  of  our  model  is  also  complicated  by  the 
existence  of  simultaneity  and  the  need,  therefore,  to  use  two-stage  least 
squares  (2SLS)  in  obtaining  final  values  for  coefficients.   We  need  a 
method,  then,  that  will: 

a)  account  for  cross-sectional  heterogeneity 

b)  correct  for  possible  autocorrelation  in  time 

c)  allow  us  to  use  2SLS,  since  the  equations  of  the  model  are 
simultaneous. 

B.3  A  Hybrid  Approach 

In  deciding  on  a  procedure  to  estimate  the  coefficients  of  this  model, 
it  must  be  kept  in  mind  that: 

a)  A  large  number  of  regressions  may  be  run  for  each  equation  in 
the  model,  using  different  specifications  and  different  regional 
breakdowns,  and  therefore  the  estimation  procedure  must  be 
convenient  and  reasonably  inexpensive, 

b)  One  purpose  of  this  model  is  to  provide  insight  into  the  structure 
and  functioning  of  the  natural  gas  industry.   Therefore,  if  the 
cross-sections  are  heterogeneous  (and  in  most  cases  we  can  expect 
that  they  will  be)  it  is  desirable  to  introduce  regional  variables 
(geological  or  economic)  to  explain  those  differences.   If  this 

is  done  properly,  most  of  the  autocorrelation  across  districts  can 
be  removed.   An  equation  for  new  discoveries  such  as  (27),  for 
example,  should  be  re-specified  in  the  form: 
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ND.  ^  =  e„  +  B.W.  ^  +  3^a.  +  e.  ^  (40) 

J,t    0    1  J,t    2  J    j,t 

where  a  is  a  geographical  "output"  variable  defined  as  in 

equation  (28) . 

c)  We  can  expect  that  in  some  equations  the  errors  will  be  autocor- 
related  in  time.   This  should  be  accounted  for. 

d)  We  are  estimating  a  simultaneous  equation  model,  and  therefore 
2SLS  will  have  to  be  used  to  ensure  that  the  estimates  are 
consistent  and  unbiased. 

Our  approach,  then,  will  be  to  use  regional  variables  whenever  possible  -- 
as  in  equation  (40)  —  to  account  for  cross-sectional  heterogeneity.   Lagged 
dependent  variables  will  be  used  to  pick  up  remaining  heterogeneity  but  will 
also  serve  to  modify  the  dynamic  response  of  the  dependent  variable.   In 
addition,  when  errors  are  autocorrelated  in  time,  we  will  use  as  a  time- 
dependent  correction  a  simple  adaptation  of  the  Hildreth-Lu  procedure.   An 
equation  such  as  (26)  would  be  tranformed  as  follows: 

^J.t  -  P^J,t-l  =  (^J,t  -  P^j.t-l>e  ^  ^j.t  -  P^j.t-1  •  <^^> 
The  value  of  the  parameter  p,  of  course,  is  not  known.  A  grid  search  can 
be  performed  on  p,  however,  and  that  value  chosen  which  minimizes  the  sum 
of  squared  residuals. 

Equation  (41)  must  in  the  end  still  be  estimated  by  two-stage  least 
squares.   The  direct  application  of  2SLS  would  lead  to  inconsistent  esti- 
mates, however,  since  the  equation  is  an  autoregressive  transformation  and 
may  also  contain  a  lagged  dependent  variable.   In  order  to  obtain  consistent 
estimates,  then,  we  will  use  Fair's  method  [7  ].   The  first-stage  regression 
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of  2SLS  is  done  as  usual,  resulting  in  a  "constructed"  series,  X    .   In 
the  second  stage,  instead  of  replacing  both  X.    and  X.   _.  with  the  con- 
structed series,  only  X    is  replaced,  so  that  OLS  is  performed  on  the 
equation: 

^j,t-P^j,t-l=  (^J,t-  ^^J,t-l)e       •  (^2) 

As  an  example,  equation  (40)  would  be  estimated  finally  by  applying  OLS  to 

^j.t  "  P^j,t-1  =  ^o^^-P^  +  "^l^^j.t  "  P"j,t-1^  "^  B2aj(l-P)  (43) 
Actually,  some  10  or  20  OLS  regressions  are  performed  on  (43),  each  of  which 
uses  a  different  value  of  p. 


4.    Estimation  of  the  Model 

One  of  the  difficulties  in  constructing  a  model  of  this  sort  is  that 
one  must  work  under  the  constraints  imposed  by  data  limitations.   Data  for 
many  variables  is  either  difficult  or  else  impossible  to  obtain,  particularly 
for  years  prior  to  1966.   In  addition,  much  of  the  data  is  extremely  noisy. 
As  a  result,  a  good  deal  of  compromise  was  often  required  in  estimating 
equations  between  functional  forms  that  are  theoretically  pleasing  and  those 
that  lend  themselves  to  the  existing  data.   This  should  be  kept  in  mind  when 
interpreting  the  estimation  results. 

A.   List  of  Variables 

A  good  deal  of  effort  was  involved  in  developing  the  data  base 

represented  by  the  variables  listed  below.   The  list,  which  includes  the 

sources  of  data,  is  divided  up  into  field  market  variables  and  wholesale 
market  variables. 
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A.l   Field  Market  Variables 

Wells 

WXO:  Exploratory  wells  completed  fcr  oil.   Consists  of  new  field  wildcats 4 

new  pool  wildcats,  outposts,  shallower  pool  tests  and  deeper  pool  tests 
(each  completed  for  oil),   oil  wells  are  defined  by  a  minimum  oil-to- 
gas production  ratio  of  (L  Ibl  oil/100,000  c.f.  gas).   Annual  data  in 
the  most  disaggregated  fonn  (cur  32  production  sub-districts)  available 
only  from  1966.   F.P.C.  Bur-au  of  Natural  Gas:   Well  Drilling  Statistics. 
Docket  No.  AR  67-1,  and  AAPG  Bulletin,  annual  wells  issue  (June  before 
1971,  July  after). 

WXG:  Exploratory  wells  completed  for  gas  or  condensate.   See  WXO  for  details. 

XF:   Total  footage  drilled  in  exploratory  wells.   Includes  gas,  oil,  conden- 
sate, and  dry.   Available  only  from  1966.   AAPG  Bulletin  1967-72  (data 
for  1966-1971). 

FWGXF:   A  geographical  variable  designed  to  act  as  a  surrogate  for  average 

drilling  costs  in  each  production  district.   It  is  the  average  explora- 
tory footage  per  gas  well  completed. 

1971    XF   . 
(FWGXF)  =j-       Z  (^r      h 

J    ^  t=1966  ^^t,j 

FWXF:  Similar  to  FWGXF,  but  this  is  the  average  exploratory  footage  per 
total  (gas  and  oil)  wells  completed.  ("Completed"  implies  that  the 
well  is  successful.) 

SR:   A  geographical  variable  defined  as  the  average  (taken  over  the  years 

1966-1971)  success  ratio  of  oil  and  gas  wells  drilled.   Defined  for 

each  district  separately. 

_  _1    ■^^'^^      ..successful  gas  wells  +  successful  oil  wells, 

3        6    ir><:<:  total  gas  and  oil  wells  drilled      .  . 

t=iyDo  J ,  t 

MMl:  Another  geographical  variable  which  is  a  surrogate  for  average  size  of 
discoveries  in  the  particular  district. 
,   1971     NU   . 

J    ^  t=1967     t-l,j 

Reserves:   All  data  from  AGA/API/CPA  Reserves  of  Crude  Oil,  etc.   Only 

available  in  disaggregated  (i.e.,  ass-diss,  non-ass.  by  production  district) 

from  1966  except  for  year-end  reserves  (YN,  YA,  US)  which  we  have  from  1965. 

In  millions  of  cubic  feet;  see  Figure  1. 


No  data  on  Louisiana  Offshore. 
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YN:   Year  end  non-associated  reserves.   Reserves  as  defined  by  the  AGA. 

YA:   Year  end  associated-dissoLv>^d  reserves.   See  YN. 

YO:   Year  end  oil  reserves. 

US:   Year  end  reserves  in  underground  storage  other  than  in  their  original 
locations. 

XN:   Extensions  of  non-associated  gas.   Includes  any  newly  proved  reserves 
already  established  in  pools  and  fields. 

XA:   Extensions  of  associated-dissolved  gas.   See  XN. 

RN:   Revisions  of  non-associated  gas.   Includes  any  proved  decreases  in  size 
of  proved  reserves  discovered  by  drilling  of  extension  wells  or  changes 
(+  or  -)  resulting  from  better  engineering  estimates  of  economically 
recoverable  reserves  in  established  pools. 

RA:   Revisions  of  associated-dissolved  gas.   See  RN. 

FN:   New  field  discoveries  of  non-associated  gas  (discovered  by  new  field 
wildcats)  . 

FA:  New  field  discoveries  of  associated-dissolved  gas.   See  FN. 

PN:  New  pool  discoveries  of  non-associated  gas. 

PA:  New  pool  discoveries  of  associated-dissolved  gas. 

XRN:  =  XN  +  RN.   Total  extensions  and  revisions,  non-associated. 

XRA:  =  XA  +  RA.   Total  extensions  and  revisions,  associated-dissolved. 

DN:  =  FN  +  PN.   Total  new  discoveries,  non-associated. 

DA:  =  FA  +  PA.   Total  new  discoveries,  associated-dissolved. 

Production:   All  data  from  AGA/API/CPA,  Reserves  of  Crude  Oil,  Etc. 

Available  in  disaggregated  form  from  1965  for  gas  (and  about  1940  for  oil). 

In  10  cubic  feet  and  10  barrels. 

QN:   Production  of  non-associated  gas.   This  is  net  production.   See  Figure  1. 

QA:   Production  of  associated-dissolved  gas.   See  QN. 

QO:   Production  of  oil. 

L:    Losses  and  waste.   API  Petroleum  Yearbook.   By  states,  Texas,  Louisiana, 
New  Mexico  undivided. 
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PG:   New  contract  price  of  interstate  sales  of  gas  at  the  wellhead.   Table  F, 
F.P.C.,  Sale  of  Natural  Gas.  Ftc. ,  1966-1970  (Cen;s/ thousand  cubic  feet). 

PO:   Wellhead  price  of  oil  ($/bbi) .   Bureau  of  Mines,  Minerals  Yearbook, 
1954-1971.   Texas  divided  into  4  sub-districts.   Otherwise,  the  same 
as  our  breakdown. 

A. 2  Wholesale  Market  Variables 

The  wholesale  market  is  presently  divided  into  two  parts,  mainline 

sales  and  sales  for  resale.   Later,  intra-state  wholesale  demand  may  be 

added.   All  variables  can  be  written  in  the  forms  A.B.C  or  B.C   .   The  first 

form  is  defined  in  terms  of  the  second  as   A.B.C  =  B.C  -  (1-A)*B.C (-1)  and 

hence  is  the  depreciated  change  in  B.C,  where  A  is  the  depreciation  rate 

(see  p.  19).   B  is  the  variable  name  (discussed  below)  and  C  refers  to  the 

number  corresponding  to  a  particular  disaggregation  of  the  United  States 

into  wholesale  markets  (see  below) . 

Variable  Names  (B) 

QI:   Mainline  industrial  sales  volume  by  interstate  pipeline  companies  by 
state  and  year  (1962-1970)  in  Mcf  (10^  cu.ft.).   From  National  Coal 
Association,  Main  Line  Natural  Gas  Sales,  1964  to  1971. 

PI:   Price  of  mainline  industrial  sales  made  on  a  firm  basis,  by  state  and 

year  (1964-1970).   Estimated  for  1962-1963  and  for  1964-1965  for  states 
with  largely  unspecified  sales.   Estimates  based  on  linear  extrapolation 
of  following  years'  data.   See  QI  for  source. 

VAM:  Value  added  in  manufacturing  in  millions  of  dollars  by  state  and  year, 
1958-1969  (1968  is  missing  and  is  interpolated).  From  U.S.  Department 
of  Commerce,  Bureau  of  the  Census,  Annual  Survey  of  Manufacturers. 

K:    New  capital  expenditures  in  millions  of  dollars.   See  VAM  for  source. 

PALT:  Price  of  alternate  fuels.   Weighted  average  (over  kilowatt-hours  generated) 
of  prices  of  fuel  oil  and  coal  consumed  by  electric  utility  industry  in 
generating  electrical  power,  by  state  and  year  (1962-1969).   Weighting 
factors  extrapolated  for  1969  from  previous  data  (1961-1968) .   From 
Edison  Electric  Institute,  Statistical  Annual  of  the  Electric  Utility 
Industry. 


In  the  sales  for  resale  equations,  we  have  converted  D.SS.6  and  D.SS.7 
to  MMcf  by  dividing  this  by  1000. 

No  data  on  Louisiana  Offshore. 
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QSR:  Sales  for  resale  volume  of  natural  gas  by  class  A  and  B  interstate 
pipeline  companies,  by  state  and  year,  1956-1970  in  Mcf.   Some  com- 
panies have  not  reported  for  1966-1968  and  some  for  1970-1971.   Data 
needs  to  be  cleaned  up.   Coupiled  from  F.P.C.  Fonu  2  reports. 

PSR:  Price  of  sale  for  resale  by  state  and  year  (see  QSR)  in  $/Mcf. 

VCC:  Value  of  construction  contracts  by  states  in  which  work  was  done, 

1956-1970,  in  millions  of  dollars.  From  Statistical  Abstract  of  the 
U.S.  and  F.W.  Dodge  Corporation,  Dodge  Construction  Contract  Statis- 
tics Service. 

Y:  Personal  income  by  state,  1956-1969,  in  millions  of  dollars.  Source: 
same  as  PSR,  and  U.S.  Bureau  of  Economic  Analysis,  Survey  of  Current 
Business. 

N:    Population  by  state  and  year,  1955-1970,  in  millions.   Source:  same 
as  above,  and  U.S.  Department  of  Commerce,  Bureau  of  the  Census, 
Current  Population  Reports. 

Wholesale  Market  Breakdowns  (C) 

C=3:  Breakdown  by  state.   Includes  twenty- two  of  the  most  important  states 
in  terms  of  volume  of  mainline  sales:   Alabama,  Arizona,  Arkansas, 
Colorado,  Florida,  Georgia,  Illinois,  Iowa,  Kansas,  Kentucky,  Louisiana, 
Minnesota,  Mississippi,  Missouri,  Nebraska,  Ohio,  Oklahoma,  Pennsylvania, 
Tennessee,  Texas,  West  Virginia,  and  Wyoming.   Others  were  omitted 
because  of  little  or  no  volume,  or  as  in  the  case  of  small  but  substantial 
states  such  as  New  Mexico,  Nevada,  Indiana  and  Michigan,  they  could  not 
be  integrated  into  a  nearby  market  because  of  extreme  price  differences. 

C=4:  Here  we  integrated  small  volume  states  in  C=3  if  they  were  contiguous 
and  had  approximately  the  same  price.   Our  breakdown  is  then 

Alabama,  Georgia 

Colorado,  Wyoming 

Kentucky,  Ohio,  Pennsylvania,  West  Virginia 

Iowa, Nebraska 

Texas ,  Oklahoma 

the  rest  of  the  states  in  C=3  being  by  themselves  market  regions. 
This  leaves  us  15  regions. 

C=6:  This  is  the  41-state  breakdown  of  sale  for  resale  which  leaves  out  only 
Alaska,  Hawaii,  Montana,  Maine,  Vermont,  North  and  South  Dakota,  and 
North  and  South  Carolina,  for  lack  of  data  (or  sales). 

C=7:  Tills  combines  several  states  in  C=6  to  get  30  districts;  i.e., 

Arkansas,  Missouri 

New  Hampshire,  Massachusetts,  Connecticut,  Rhode  Island 

New  Mexico,  Utah,  Colorado 

Nevada,  Oregon,  Washington 

Virginia,  Maryland 

Wyoming,  Idaho. 
Delaware  is  deleted  because  its  price  is  too  diff-arent  from  other  states 
nearby  to  justify  including  its  small  volume  of  sales. 


B.   Statistical  Results 

All  of  the  regression  results  described  below  were  obtained  using 
ordinary  least  squares  (OLS) .   It  was  not  necessary  to  use  two-stage  least 
squares  (2SLS)  —  unless  otherwise  indicated  —  because  the  final  equation 
forms  that  best  fit  the  data  contained  no  endogenous  variables  that  were 
not  lagged.   This  was  not  the  result  of  a  predisposition  to  eliminate  simul- 
taneity from  the  model,  but  merely  a  statistical  result. 

In  addition,  we  did  not  use  tie  modified  version  of  Hildreth-Lu  des- 
cribed in  section  3.B  in  obtaining  these  results.   It  was  our  feeling  that 
this  initial  effort  at  estimating  the  model  should  be  preliminary  and  ex- 
perimental, designed  more  to  test  the  data  and  the  variety  of  relationships 
that  might  fit  the  data.   We  do  expect  to  use  that  method,  however,  in 
estimating  a  future  version  of  this  model. 

All  of  our  equations  are  linear  in  form  with  the  exception  of  the 
wells  equation  (for  WXG)  which  is  in  log-log  form.   The  log  version  of 
this  equation  fit  far  better  than  did  any  of  a  whole  variety  of  linear 
versions.   It  is  interesting  to  note  that  for  all  of  the  other  equations 
of  the  model  none  of  the  log  versions  fit  as  well  as  the  linear  versions. 
Future  work  on  this  model  will  include  the  testing  of  other  equation  forms, 
including  some  that  are  non-linear.   Also,  we  expect  to  be  able  to  obtain 
reserves  data  for  the  years  1964  and  1965  (our  present  reserves  data  only 
covers  the  years  1966-1971),  so  that  the  supply  equations  for  the  field 
market  can  be  estimated  over  a  longer  time  period. 
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B.l  Field  Market  Equations 

The  field  market  portion  of  our  model  contains  eight  stochastic 
equations  together  with  an  identity  to  define  changes  in  reserves.   Seven 
of  the  stochastic  equations  (wells,  non-associated  discoveries,  associated 
discoveries,  non-associated  extensions,  associated  extensions,  non-associated 
revisions,  and  association  revisions)  determine  the  supply  of  new  reserves, 
while  the  last  (explaining  the  price  level  before  1962)  determines  demand. 

The  price  equation,  of  course,  can  only  be  estimated  using  data  before 
1962.   We  were  able  to  obtain  data  for  the  price  equation  beginning  in  1956, 
and  the  cross-sections  were  very  heterogeneous.   As  a  result,  our  initial 
estimation  results  for  this  equation  were  rather  unsatisfactory,  and  rather 
than  extrapolate  a  poorly  fitted  equation  from  1962  to  1972  (and  then  to 
1980  for  a  forecast),  we  decided  to  leave  the   equation  out  of  this  preli- 
minary version  of  the  model.   The  wellhead  price  is  exogenous  anjrway  after 
1962,  and  therefore  leaving  this  equation  out  does  not  affect  our  simulation 
experiments  and  forecasts  for  new  reserves  and  production  in  the  future.   It 
does  make  it  impossible,  however,  to  obtain  a  direct  estimate  of  excess  de- 
mand in  the  field  market.   An  indirect  estimate  can  be  obtained  though  by 
equilibrating  the  wholesale  market  with  the  field  market  under  the  assump- 
tion of  some  minimum  required  reserves-to-production  ratio.   This  point  will 
be  discussed  in  more  detail  later. 

Our  final  equation  for  wells  is  in  log  form  and  relates  completed  ex- 
ploratory wells  to  the  wellhead  price  of  gas  (which  is  the  main  incentive  to 
drill),  the  production  of  oil  (also  an  incentive  to  drill,  since  oil  may  be 
discovered),  and  previous  drilling.   The  equation  is  shown  below  together  with 
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2  , 

the  R  ,  F-statistic,  standard  error  of  the  regression,  Durbin-Watson  statistic, 

and  statistics  (in  parentheses). 

log  WXG  =  -3.335  +  1.022  PG  _  +  0.1212  QO  _  +  0.6684  WXG  _  (44) 

(-2.10)   (2.31)    ^   ^      (1.61)    ^  (6.21)     ^   ^ 

R^  =  .7096         F  =  40.7        S.E.  =  0.591       DW  =  1.75 
Other  versions  of  this  equation  were  estimated  both  in  log  form  and  in 
linear  form  (none  of  the  linear  forms  were  satisfactory) ,  and  the  results 
are  summarized  in  Table  4.   It  is  interesting  to  note  that  the  geographical 
variable  MMl  which  measures  the  "richness"  of  each  production  district  was 
not  significant  in  any  of  the  versions,  and  as  a  result  was  dropped  from  the 
final  form.   Other  specifications  not  listed  in  Table  4  were  tested  but  were 
found  to  be  inferior.   It  was  thought  that  lagged  new  discoveries  could  he 
used  as  an  explanatory  variable  sinte  recent  discoveries  in  a  region  might 
stimulate  further  drilling,  but  the  relationship  was  not  found  to  be  signi- 
ficant. 

Our  equation  for  new  discoveries  of  non-associated  gas  (DN)  does  contain 
an  unlagged  endogenous  variable,  namely  wells,  WXG,  and  therefore  was  esti- 
mated using  two  stage  least  squares  (a  fitted  series  for  WXG  was  formed  by 
regressing  it  on  a  set  of  instrumental  variables).   Our  geographical  variable 
MMl  (the  average  past  size  of  discoveries  per  well  on  a  regional  basis)  was 
used  in  this  equation  and  proved  quite  significant.   This  final  form  of  the 
2SLS  regression  was: 

DN  =  -141773.-  +  25.085  MMl  +  4189.3  WXG  +  0.261  DN  (45) 

(-2.66)     (3.11)        (1.78)      (1.54) 

R^  =  .721      F  =  43.0       S.E.  =  2.2  x  lo^    DW  =  1.16 
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Table  4:   Alternative  Wells  Equations 
(t-statistics  in  parentheses) 


A. 

LOG  FORM 

r2 

const. 

PG 

^^t-1 

QO 

^°t-l 

^°t-l 

MMl 

^^-1 

.695 

-1.693 
(-1.11) 

0.6055 
(1.10) 

0.4568 
(0.340) 

0.0111 
(0.129) 

0.7640 
(7.719) 

.711 

-3.096 
(-1.86) 

1.015 
(2.27) 

0.1367 
(1.69) 

-.047 
(-.547) 

0.6734 
(6.19) 

.711 

-3.063 
(-1.85) 

1.019 
(2.267) 

0.132 
(1.66) 

-.046 
(-.532) 

0.6706 
(6.0,2) 

.709 

-3.298 
(-2.08) 

1.025 
(2.29) 

0.118 
(1.584) 

0.666 
(6.093) 

.693 

-2.690 
(-1.44) 

0.8265 
(1.54) 

0.103 
(1.312) 

0.6905 
(5.854) 

LINEAR  FORM 


PG 


t-1 


QO 


QO 


t-1 


PO 


MMl 


t-1 


.552 

-9.43 
(-0.85) 

0, 
(1. 

.985 
.61) 

.0000168 
(.855) 

.000573 
(1.26) 

0.388 
(2.81) 

.545 

-5.812 
(-0.22) 

0.968 
(1.08) 

-0.910 
(-.08) 

.000765 
(1.91) 

0.436 
(3.31) 

.554 

-8.14 
(-0.79) 

0.942 
(1.60) 

.000021 
(0.99) 

.000546 
(1.20) 

0.393 
(2.90) 
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Other  forms  of  this  equation  use  a  first  difference  in  wells,  multiply  wells 
by  the  geographical  variable  MMl,  and  include  the  price  of  gas,  but  these 
did  not  fit  as  well.   An  attempt  was  also  made  initially  to  divide  DN  up  into 
pools  (PN)  and  fields  (FN) ,  but  we  were  unable  to  obtain  an  equation  that  had 
a  reasonable  fit  for  the  latter. 

Our  equation  for  associated  discoveries  (DA)  also  contained  the  wells 
variable  WXG  unlagged,  and  was  therefore  also  estimated  using  two  stage  least 
squares.   The  geographical  variable  MMl  was  not  significant  in  this  equation, 
but  that  was  not  surprising  since  the  "payoff"  for  associated  discoveries  is 
linked  so  closely  to  both  oil  discoveries  and  oil  production  (many  oil  wells 
also  contain  gas  which  can  be  tapped  after  much  of  the  oil  has  been  produced) . 
As  a  result,  we  included  as  an  exogenous  variable  in  this  equation  the  pro- 
duction of  oil  QO.   The  final  form  of  the  equation  is: 

DA  =  -1078A.2  +  0.0685  QO  +  381.64  WXG  +  0.426  DA  (46) 

(-2.78)    (3.40)        (2.33)      (4.01) 

R^  =  .833      F  =  83.2       S.E.  =  1.5  «  10^    DW  =  1.73 

Again  an  attempt  was  made  to  divide  associated  discoveries  into  those 
from  pools  (PA)  and  those  from  fields  (FA) : 

PA  =  -4984.8  +  .0436  QO  +  76.35  WXG    +  0.586  PA  (47) 

(-1.68)    (0.66)      (0.71)        (6.14) 

R^  =  .843      F  =  89.7       S.E.  =  1.25  ^  lO'^   DW  =  1.94 

FA  =  1002.3  +  .0253  QO  +  18.90  WXG  _  +  0.104  FA  (48) 

(0.668)   (3.35)      (0.36)    ^   ^      (0.74)    '^"■ 

R^  =  .368      F  =  9.71       S.E.  =  7.0  ^  10"^    DW  =  1.80 
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As  can  be  seen,  these  results  are  not  as  satisfactory,  much  of  the  difficulty 
arising  '^rom  the  instability  in  the  series  for  FA.   Thus,  the  equation  in 
the  aggregated  form  (i.e.,  for  DA)  was  used. 

Extensions  of  non-associated  gas  (XN)  are  related  positively  to  the 
price  of  gas,  negatively  to  the  price  of  oil  (since  a  higher  oil  price  is 
a  disincentive  for  companies  which  explore  for  both  oil  and  gas  to  spend 
money  extending  gas  fields),  and  positively  to  the  number  of  wells  drilled: 

XN  =  8.816  X  10^  +  57813.0  PG  -  6.1  ^  10^  PO  +  5269.8  WXG    +  0.615  XN  _ 
(1.16)         (2.18)      (-1.81)  (1.99)         (10.96)  ^ 

R^   =  .774      F  =  42.0       S.E.  4.4  x  lo^      DW  =  1.54       ^^^■' 

Note  that  it  is  not  necessary  to  use  2SLS  to  estimate  this  equation  since  the 

price  of  gas  PG  is  exogenous  after  1962,  i.e.,  during  the  period  over  which 

the  equation  is  estimated. 

Extensions  of  associated  gas  (XA)  depend  on  the  production  of  oil  QO 

(dissolved-associated  gas  fields  are  extended  to  allow  for  increased  oil 

production),  and  also  the  previous  year's  new  discoveries  of  associated  gas: 

XA  =  24521.9  +  0.4195  QO  +  0.2671  DA    +  0.2327  XA  (50) 

(-2.07)   (4.36)      (0.60)         (1.62) 

R^  =  .757      F  =  52.0       S.E.  =  6.3  ^  lo'^    DW  =  1.65 

Revisions  of  non-associated  gas  (RN)  depend  positively  on  the  previous 
year's  total  reserves  of  non-associated  gas  (YN) ,  and  negatively  on  the 
change  in  the  previous  year's  reserves  of  non-associated  gas.   Basically, 
large  short-run  increases  in  reseirves  usually  limit  revisions  in  the  fol- 
lowing year  (which  apply  more  to  established  discoveries),  but  in  the  long 
run  a  higher  level  of  reserves  results  in  a  higher  average  level  of  revisions. 
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Our  equation  is  : 

RN  =  -1.946  X  10^  +  0.01A3  YN  _  -  0.1726  AYN  _  +  0.4946  RN  (51) 

(-1.71)        (2.11)      ^  (-4.09)     ^   ^      (3.24)    ^  ^ 

R^  =  .511      F  =  17.4       S.E.  =  6.03  ^  lo^   DW  =  1.52 

Note  that  the  coefficient  of  the  autoregressive  term  (RN  _^ )  is  close  to 
.5,  and  thus  the  long  run  impact  on  revisions  resulting  from  an  increase 
in  reserves  of  1  Mcf  is  approximately  equal  to 
1 


1  -  .5 


(.0143)  =  .0286  Mcf 


increase  in  the  stock  of  non-associated  reserves  implies  about  a  .03  Mcf 
per  year  increase  in  the  flow  of  non-associated  revisions. 

Associated  revisions  (RA)  are  extremely  erratic  and  difficult  to  ex- 
plain using  a  simple  linear  regression  model.   There  seems  to  be  no  rela- 
tionship between  this  variable  and  lagged  reserves  of  associated  gas,  mostly 
because  associated  revisions  are  linked  more  closely  to  oil  reserves  (which 
we  do  not  include  in  this  model) .   We  modeled  associated  revisions  by  re- 
lating it  to  the  previous  year's  associated  discoveries  and  the  previous 
year's  associated  extensions,  although  the  relationship  is  somewhat  dubious: 

RA  =  -99078.6  +  4.766  DA  _  +  0.844  XA  _  -  0.0543  RA  _  (52) 

(-1.30)    (1.56)    ^  ^      (0.89)    ^   ^    (-0.39)    ^  ^ 

R^  =  .236      F  =  5.2        S.E.  =  4.73  x  10^   DW  =  1.47 


Note  that  AYN  _   =  YN  _   -  YN  _  . 
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B.2  Wholesale  Market  Equations 

There  are  currently  four  stochastic  equations  in  our  representation 
of  the  wholesale  market  for  gas.   Consumption  demand  equations  are  estimated 
for  mainline  sales  (QI)  and  for  sales  for  resale  (QSR) ,  and  price  markup 
equations  are  also  estimated  for  both  mainline  sales  (PI)  and  sales  for  re- 
sale (PSR) .   Note  that  we  are  describing  the  market  only  for  interstate  sales 
of  gas,  and  the  production  identity 

Q  =  QI  +  QSR  (53) 

defines  only  interstate  production.   This  must  be  kept  in  mind  when  Q  is  used 
later  to  calculate  the  reserves-to-production  ratio. 

Before  the  mainline  demand  equation  could  be  estimated,  an  estimate  had 
to  be  made  of  the  depreciation  rate  for  industrial  gas-burning  equipment  (r) . 
This  was  done  by  estimating  the  rate  simultaneously  with  the  other  coefficients 
of  a  total  industrial  demand  equation  that  did  not  contain  an  autoregressive 
term: 

QI  =  (1-  0.104)  QI    +  17175.8  -  97554.8  PI    +  9144.6  FIR  (54) 

(2.95)  (0.86)     (-1.62)         (1.18) 

+  1.576  K    +  79150.3  PALT 
(0.46)        (2.88) 

R^  =  .317      F  =  7.8  DW  =  2.11 

Here  FIR  is  the  f irm-to-interruptible  ratio,  K  is  new  capital  expenditures, 
and  PALT  is  a  price  index  of  alternative  fuels.   Note  that  the  depreciation 
rate  r  came  out  to  be  about  10%  and  was  significant  (t  =  2.95).   A  deprecia- 
tion rate  of  0.1  was  then  used  to  calculate  our  "new"  demand  variable 
6QI  =  QI  -  (l-r)QI  _   for  use  in  the  final  equation  for  mainline  demand. 
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Our  mainline  demand  equation  contains  the  same  explanatory  variables 
as  equation  (54)  above,  except  that  it  also  contains  an  autoregressive  term. 
This  is  added  to  help  explain  the  dynamic  response  of  new  demand  (i.e.,  it 
imposes  a  Koyck  distribution  on  the  explanatory  variables)  and  also  to  ac- 
count partially  for  whatever  heterogeneity  might  remain  in  the  cross- 
sections.   Our  final  equation  was: 

6QI  =  10951.0  -  74609.8  PI    +  7611.9  FIR  _  +  1.063  K  _  (55) 

(0.64)     (-1.69)     ^        (1.10)    ^  ^      (0.35)  ^   ^ 

+  69467.0  PALT    +  0.181  6QI 
(2.85)        ^      (1.62)    " 

R^  =  .325      F  =  8.1        S.E.  =  1.8  ^  lo^    DW  =  2.11 

2 
Note  that  the  R  and  t-statistics  are  low,  but  this  reflects  the  scaling  of 

the  dependent  variable.   If  the  depreciation  rate  r  was  set  to  0.18  instead 

2 
of  0.10,  the  mean  of  the  dependent  variable  would  be  larger  and  the  R   also 

larger.   The  regression  result  using  r  =  0.18  is  shown  below: 

6QI  =  24418.3  -  1.30  x  lO^  PI    +  13631.6  FIR    +  3.154  K  _         (56) 
(1.38)   (-2.73)        ^  ^        (2.05)  (1.01)  ^   ^ 

+  85583.6  PALT    +  0.290  6QI 
(3.40)        ^   (2.73)    ^   ^ 

R^  =  .569      F  =  22.2       S.E.  =  1.8  x  lo^    DW  =  2.14 

4 
Note  that  the  standard  error  in  both  cases  is  the  same  (1.8  x  10  ) . 

An  attempt  to  estimate  the  depreciation  rate  for  the  sales  for  resale 

demand  equation  was  unsuccessful,  since  the  rate  that  resulted  was  negative. 

We  therefore  used  the  same  depreciation  rate  (0.10)  as  was  used  in  the  main^ 

line  demand  equation.   Explanatory  variables  include  the  wholesale  price. 
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the  price  of  oil  (oil  is  the  main  competing  fuel  for  residential  and  commer- 
cial demand)  and  total  personal  income  (on  a  state-wide  basis).   The  autore- 
gressive  term  was  also  included.   The  time-bounds  for  this  equation  went 
from  1962  to  1969,  and  30  districts  were  used  (whereas  an  aggregated  set  of 
15  districts  was  used  in  the  mainline  demand  equation  —  see  Table  4).   The 
resulting  equation  is  shown  below: 

6QSR  =  6.32  X  10^  -  1.762  ^  10^  PSR  +  1.121  x  lo^  PO  +  15599.5  Y      (57) 
(0.3C)      (-3.67)  (4.95)  (8.07) 

+  0.103  6QSR 
(1.46) 

R^  =  .312      F  =  26.6       S.E.  =  6.35  ^  lO''   DW  =  2.38 

2 
Again,  the  low  R  reflects  the  dimensioning  of  the  dependent  variable. 

In  keeping  with  the  arguments  on  page  23,  the  wholesale  price  markup  was 
seen  to  depend  on  volumetric  capacity  and  average  mileage  from  production 
regions  to  the  particular  wholesale  market  region.   Both  of  these  explanatory 
variables  are  assumed  to  be  roughly  constant  in  time,  so  that  our  price  mark- 
up regressions  are  actually  cross-section  regressions  resulting  in  constant 
markups.   For  mainline  sales,  the  markup  is  given  by: 


(58) 


(3.71)   (23.46)         (-3.25) 

R^  =  .861      F  =  463.6      S.E.  =  2.71         DW  =  0.17 

Note  that  the  coefficient  of  volumetric  capacity  v  is  negative,  since  a 

larger  capacity  implies  lower  average  costs.   For  sales  for  resale,  our 

markup  equation  is  given  by: 

PSR  =  PG  +  13.672  +  0.5973  M  -  4.828  x  lo"'^  v 
(14.82)   (13.88)    (-13.56) 

R^  =  .688      F  =  426.9      S.E.  =  5.05         DW  =  0.23 


(59) 
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This  completes  the  estimation  of  the  model.   In  the  next  section, 
the  equations  are  combined  (together  with  the  appropriate  accounting  iden- 
tities), and  the  complete  model  is  simulated. 

5.    Use  of  the  Model  for  Policy  Simulations 

Before  the  model  can  be  simulated  into  the  future,  a  mechanism  must  be 
specified  by  which  gas  is  distributed  from  production  regions  to  consumption 
regions,  since  the  equations  as  they  now  stand  cannot  be  simply  combined  and 
simulated  directly.   The  wellhead  price  that  appears  in  the  field  market 
equations,  for  example,  applies  to  a  different  set  of  regions  than  does  the 
wellhead  price  that  appears  in  our  wholesale  price  markup  equations.  Similarly, 
additions  to  reserves  (from  discoveries,  etc.)  are  made  in  field  market  re- 
gions, while  depletions  of  reserves  (resulting  from  production)  are  indexed 
by  wholesale  regions.   Furthermore,  the  wholesale  regions  themselves  are 
different  for  mainline  sales  and  sales  for  resale.   Thus  the  model  cannot 
be  equilibrated  unless  there  is  some  means  by  which  reserve   additions  from 
field  markets  can  be  added  together  and  then  redistributed  to  wholesale  markets. 

To  do  this  we  have  divided  the  U.S.  into  three  "zones"  and  then  postulated 
three  large  regional  field  markets  (East,  North,  and  West),  in  which  gas  pro- 
duced in  the  various  production  districts  can  be  accumulated  and  then  distri- 
buted to  the  various  regional  markets  for  mainline  sales  and  sales  for  resale. 
The  regional  breakdown  that  we  have  selected  is  shown  in  Table  5.   Note, 
for  example,  that  gas  from  the  Texas  Gulf  Coast,  East  Texas,  and  Louisiana 
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Table  5:   Regional  Interrelationships  for  Simulation 


Regional 

Mainline     i 

Sales  for 

Production  Districts 

Field  Markets    | 

Markets 

Resale  Markets 

Texas  RRC  //I      "^ 

Alabama, 

Alabama 

Texas  RRC  #2 

Georgia 
Florida 

New  England 
Florida 

Texas  RRC  //3 

Georgia 

Texas  RRC  //A 

Appalachian 

Kentucky 

Texas  RRC  //5,6     > 

East 

Louisiana 

Louisiana 
Maryland,  Virginia 

Louisiana  North 

Mississippi 

Mississippi 

Louisiana  South 

New  Jersey 

(including  Off- 
shore)        J 

New  York 

Pennsylvania 

Tennessee 

Tennessee 

West  Virginia 



Arkansas , 
Missouri 

Illinois 
Iowa, Nebraska 

Arkansas,  Missouri 

Indiana 

Illinois 

Iowa 

Texas  RRC  //7B,8A,9  ' 

\ 

Kansas 

Kansas 
Michigan 

Texas  RRC  #10 

Minnesota 

Minnesota 

* 

North 

Nebraska 

Oklahoma 

Ohio 

Oklahoma, 

Oklahoma 

Kansas 

Texas 

Texas 

J 

1 

Wisconsin 

Texas  RRC  //7C,  8 

t 

i 

1 
1 

1  Arizona 

Arizona 

New  Mexico  Northwest  . 

1 

! 

1   West 

1  Colorado, 

1   Wyoming 

California 

Rocky  Mountain 

New  Mexico  Southwest 

1 

Idaho,  Wyoming 

/  i 

Nevada,  Oregon,  Washing 

ton 

Montana,  Colorado,  Utalj   oth  r 
Calif,  Wyoming,  West   j 

1 

VA^.-.i^;  ., : 

1 
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is  sold  to  pipelines  which  begin  in  the  Regional  Field  Market  East  and  then 
transmit  the  gas  east  to  New  England,  the  Appalachian  region,  the  Eastern 
Seaboard,  and  the  Southeast.   Gas  from  North  Texas,  Oklahoma,  and  Kansas 
is  sold  in  the  North  Regional  Field  Market  to  pipelines  distributing  gas 
to  the  North  Central  and  South  Central  U.S.   Gas  produced  in  West  Texas  and 
New  Mexico  is  sold  in  the  West  Regional  Field  Market  to  pipeline  companies 
which  distribute  it  to  the  West  and  Northwest  U.S.   Since  gas  produced  in 
Montana,  Colorado,  Utah,  California,  Wyoming,  and  West  Virginia  does  not 
contribute  significantly  to  Interstate  mainline  sales  or  sales  for  resale, 
we  have  omitted  them  from  the  interactive  part  of  the  model,  though  they 
are  still  included  in  terms  of  their  contribution  to  year-end  reserves  through 
new  discoveries,  extensions,  revisions,  etc. 

Note  that  this  breakdown  of  field  markets  need  not  correspond  to  actual 
economic  markets.   We  have  not  yet  been  able  to  statistically  identify  the 
geographical  boundaries  of  actual  field  or  wholesale  markets,  and  thus  for 
purposes  of  simulation  we  are  simply  trying  to  aggregate  districts  in  a 
"reasonable"  way.   An  extension  of  this  research  will  include  an  attempt 
at  better  identifying  regional  markets  and  their  interrelationships.   As 
we  have  explained  earlier  in  this  paper,  the  interaction  of  regional  markets 
is  an  important  aspect  of  the  natural  gas  industry. 

A.    Initial  Policy  Simulations 

In  the  policy  simulations  that  follow,  we  have  not  included 
equation  (49)  for  non-associated  extensions,  since  it  is  unrealistic  in 

explaining  the  behavior  of  this  variable  when  the  price  of  gas  at  the  wellhead 
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is  fixed  to  be  constant  and  the  price  of  oil  is  growing.   Thus,  in  these 
initial  simulations,  non-associated  extensions  (XN)  is  taken  to  be  exogenous 
and  allowed  to  grow  at  5%  from  its  1970  value. 

Table  6  shows  our  assumptions  about  the  growth  of  the  exogenous 
variables.   The  price  of  oil  is  assumed  to  grow  at  6%,  which  is  probably 
on  the  high  side  and  is  the  cause  of  the  poor  performance  of  the  XN  equation 
(with  oil  prices  rising  that  rapidly  and  gas  prices  fixed,  it  does  not  pay 
to  extend  non-associated  gas  fields).   In  the  next  set  of  simulations,  we 
will  reduce  the  growth  rate  of  the  price  of  oil  to  2%  per  year,  and  then 
include  the  XN  equation  in  the  model. 

We  ran  six  simulations  using  this  version  of  the  model  together  with 
the  exogenous  variable  forecasts  summarized  in  Table  6.   In  the  first  simu- 
lation, we  hold  the  ceiling  price  of  gas  fixed  at  its  1970  values  (on  a 
region-by-region  basis).   In  the  second,  we  allow  the  ceiling  price  to  in- 
crease by  one  cent  per  year,  so  that  in  1980  it  is  ten  cents  higher  than  its 
1970  level  in  each  production  discrict.   In  the  remaining  simulations  it 
increases  by  2c  per  year,  3c  per  year,  4c  per  year,  and  5c  per  year.   The 
results  are  summarized  graphically  for  several  key  aggregate  variables  in 
Figures  5  through  10,  and  on  an  East-West-North  regional  basis  in  Figures 
11  through  15. 

Note  that  the  drilling  of  wells  is  extremely  sensitive  to  the  ceiling 
price  of  gas,  and  therefore  additions  to  reserves  (through  new  discoveries 
and  extensions  and  revisions)  is  indirectly  sensitive,  although  with  a  time 
lag  of  about  two  years.   With  no  change  in  the  ceiling  price  of  gas  during 
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Table  6:   Forecast  Assumptions  for  Exogenous  Variables 


Price  of  Oil  (PO) 

Quantity  of  Oil  (QO) 

Price  of  Alternate  Fuels  (PALT) 

Personal  Income  by  State: 

East: 
North: 

West: 

Capital  Expenditures  (K) 
(Business  Fixed  Investment) 


grows  at  6%  per  year 

grows  at  6%  per  year 

grows  at  6%  per  year 

grows  at  7.2%  per  year 

grows  at  7.5%  per  year 

grows  at  7.7%  per  year 

grows  at  6%  per  year 


Firm-to-Interruptible  Ratio  (FIR)   assumed  constant 
Mileage  (M)  assumed  constant 

Volumetric  Capacity  (v)  assumed  constant 

Non-Residential  Structures         grows  at  3.5%  per  year 
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the  decade  1970-1980,  the  number  of  wells  drilled  declines  slightly,  total 
new  discoveries  remain  about  constant,  and  extensions  and  revisions  decline 
rapidly.   Total  demand  Jor  gas  (i.e.,  production),  however,  increases  by 
about  11%  per  year,  for  a  total  increase  of  150%  over  the  decade.   This 
increase  in  production  demand  is  to  be  expected  since  the  wholesale  price 
is  constant  (the  markup  is  constant),  while  the  price  of  alternate  fuels, 
GNP,  and  capital  expenditures  are  all  increasing.   If  this  demand  were  to 
be  somehow  satisfied  (by  depleting  reserves),  we  would  run  out  of  reserves 
by  1980.   This  demand,  of  course,  would  not  be  satisfied.   By  the  time  the 
reserves-to-production  ratio  dropped  to  about  7  or  8  (i.e.,  by  1974), 
considerable  rationing  would  occur,  and  many  consumers  would  simply  not 
be  able  to  purchase  gas. 

Even  with  a  ceiling  price  increase  of  2c  every  year,  demand  for  reserves 
would  increase  faster  than  the  supply,  and  the  reserves-to-production  ratio 
would  drop  to  about  2.3  by  1980.   Many  more  wells  would  be  drilled  in  this 
case,  and  both  new  discoveries  and  extensions  and  revisions  would  increase 
accordingly,  but  this  increase  in  reserves  would  be  outstripped  by  the  more 
rapidly  growing  increase  in  production.   Even  with  a  2c  per  year  increase  in 
the  price  of  gas,  increases  in  the  prices  of  other  fuels,  in  GNP,  and  in 
capital  expenditures  would  result  in  a  doubling  of  wholesale  gas  demand  by 
1980. 

A  ceiling  price  increase  of  4c  per  year  or  5<;  per  year  would  result  in 
increases  in  reserve  additions  that  would  be  larger  than  the  increases  in 
wholesale  demand,  so  that  the  reserves-to-production  ratio  would  increase 
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over  the  decade,  reaching,  by  1980,  23  for  a  Ac  annual  price  increase  and  41 
for  a  5C  increase.   This  is  the  tesult  of  a  dramatic  increase  in  well  dril- 
ling together  with  an  increase  ii  wholesale  demand  of  only  4%  to  6%  per 
year. 

If  F.P.C.  objectives  are  to  keep  the  reserves-to-production  ratio 
approximately  constant  at  around  12,  then  these  results  indicate  that 
the  ceiling  price  should  be  increased  by  about  3-3. 5C  per  year.   In  this 
case,  additions  to  reserves  would  be  commensurate  with  increases  in  produc- 
tion demand,  and  the  reserves-to-production  ratio  would  remain  about  the 


Note  that  all  of  these  results  make  certain  assumptions: 

(1)  That  exogenous  variables  will  grow  at  certain  rates  (for  example,  the 
price  of  oil  is  assumed  to  grow  at  6%  per  year) . 

(2)  Geological  conditions  will  remain  the  same,  i.e.,  we  will  not  simply 
begin  to  run  out  of  gas.   Thus  we  assume  that  wells  drilled  in  1975 
will  produce  as  much  gas  as  those  drilled  in  1970. 

(3)  We  will  not  run  out  of  oil  or  see  any  rationing  of  oil. 

(4)  The  basic  determinants  of  the  demand  for  gas  will  not  change  in  time. 
This  last  assumption  may  be  the  most  questionable,  at  least  in  the 

short-run.   We  may  have  begun  to  see  in  the  past  year  or  two  a  change  in 
"ecological  consciousness"  in  the  U.S.   On  the  whole,  gas  is  a  somewhat 
cleaner  fuel  than  oil,  and  this  may  result  in  additional  increases  in  the 
demand  for  gas. 
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Figure   5:      Total   Wells   Drilled   and   Completed 
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Figure  6:   Total  New  Discoveries 
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Figure  7:   Total  Extensions  and  Revisions 
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Figure  8:   Total  Year-End  Reserves 
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Figure   9:      Total   Production 
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Figure  10:   Reserves-to-Production  Ratio 


SIMULATIOtI  CHJTPUT  BY  l>^IABLE 
RP. RATIO  -  DtFINlTIO< 


45 


20 


1970 
TIt1E  BOUNDS   1970  TO  1988 


1975 


1989 


Numbers  at  far  right  indicate  the 
price  increase  in  cents  per  year. 


Fifture    11:      Wells   Drilled    and    Co~iplctcd  ,    by    Region 
Wells    -    East 


69 


R    0 


1970 


1975 


1980 


Wells    -    West 


IS^ 


Wells    -   North 


4«eo 


2000 


Figure    12:      New  Discoveries,    bv   Resion 
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Figure    13:       Extensions    anJ    Revisions,    By    Rc>;lon 
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Figure    14:      Year-End   RcservfS ,    By  Region 
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Figure    15:      Production,    By   Regions 
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6.    Tentative  Conclusions 

The  process  of  construction  of  an  econometric  model  of  natural  resource 
utilization  is  iterative.   This  is  a  report  on  the  first  of  the  iterations; 
as  improvements  are  made  in  the  data  base,  by  extending  the  coverage  and  in- 
creasing the  accuracy  of  the  observations,  the  model  will  be  refitted.   As 
mistakes  in  the  formal  structure  of  the  model  are  revealed  by  our  readers, 
there  will  be  additional  computations  of  the  equations.   And  as  the  iterations 
occur,  we  expect  the  findings  to  change. 

Given  the  first-round  results  in  the  preceding  section,  we  can  formulate 
preliminary  conclusions,  however.   Following  a  strategy  of  well-publicized 
price  increases,  the  Federal  Power  Commission  can  eliminate  the  domestic  gas 
shortage  by  1980.   The  best  of  the  (feasible)  price  changes  would  be  3  cents 
per  annum.   In  contrast,  1  cent  per  annum  would  eliminate  the  entire  stock 
of  domestic  reserves,  and  5  cents  per  annum  would  not  be  maintainable  in  the 
markets  for  gas  because  there  would  be  a  veritable  flood  of  additional  dis- 
coveries relative  to  additional  production  demand.   The  3  cent  per  annum  in- 
crease would  "find"  enough  gas  while  restraining  demands  for  a  domestic  12:1 
reserve-to-production  ratio  —  close  to  that  before  the  shortage.   The 
resulting  level  of  wholesale  prices  would  be  less  than  that  which  would 
occur  under  a  continuation  of  the  present  ceilings  and  importation  of  liqui- 
fied gas  at  $1.15  per  Mcf. 

The  "public  strategy"  set  out  here  for  solving  the  gas  "crisis"  is  one 
of  gradual,  publicized  increases  in  F.P.C.  ceiling  prices.   An  alternative 
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would  be  to  deregulate,  In  stages  or  abruptly,  so  that  markets  cleared 
through  price  increases  in  the  middle  1970' s.   We  cannot  simulate  this 
alternative  strategy  at  the  present  time,  since  this  first  model  was 
built  on  1960-1970  conditions  in  which  prices  were  exogenous.   An  attempt 
will  be  made  to  anticipate  the  effects  of  endogenous  prices,  however,  in 
forthcoming  iterations.   This  will  be  done  by  formulating  and  testing  a 
functional  relation  for  "demand  for  reserves"  based  on  1955-1960  experience 
before  ceiling  prices;  the  addition  of  this  function  to  the  model  "closes" 
the  system  and  determines  prices  at  equilibrium  levels  of  reserves  and  pro- 
duction. 

This  approach  will  be  most  tentative.   Extrapolating  behavioral  relations 
fitted  on  data  from  the  1950 's  is  likely  to  lead  to  large  forecast  errors. 
But  it  may  provide  insights  into  the  contrasting  balance  of  short-term 
market  adjustments  —  contrasting  to  the  "lock-step"  price  increases  posited 
in  the  simulations  of  regulation  provided  here. 

The  case  for  less  price  control  is  made  here,  in  these  initial  finding!} 
of  gains  for  consumers  from  3  cent  per  annum  price  increases.   The  case  for 
deregulation  is  made  on  philosophical  grounds  —  we  should  not  regulate 
effectively  competitive  markets,  or  use  price  controls  to  solve  problems 
in  taxing  the  windfall  gains  of  resource  producers.   Improvements  in  this 
work  may  lead  to  quantitative  findings  which  indicate  that  immediate 
deregulation  is  preferred. 
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